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ABSTRACT 



The ability to make computer images more realistic is becoming more 
important a s the hardware for producing such images is becoming less expensive 
and hence more available. The key to producing realistic images lies in the 
algorithms that can take full advantage of the hardware to produce them. In this 
study, we look at a prototype of a ray tracer, as presented in [Ref. 1]. Ray tracing, 
in combination with a global illumination model, currently provides the most 
realistic images that can be generated on general purpose computing hardware. 
The prototype was succesfully implemented on an IBM AT clone. 
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I. INTRODUCTION 



From the beginning of recorded history mankind has always had the need to 
create pictures. The reasons for creating these pictures ranges from the aesthetic, 
pretty pictures are nice to look at, to the functional, pictures can be an excellent 
way to communicate information. As mankind progressed, so did his ability to 
create pictures, although the techniques used to create pictures basically stayed 
the same. The advent of computers gave man yet another tool with which to 
create pictures. 

The rapid increase in technology has made computer graphics a rapidly 
growing field. For the first time since man started drawing pictures, completely 
new techniques needed to be developed. In computer graphics the brush, paint 
and canvas are replaced by the mouse, algorithm and display. However even 
though the tools have changed, the same problems remain: how to make the 
picture look better, be it either more pleasing to the eye or to get the information 
across more clearly. 

Two of the most common and difficult problems in computer graphics are the 
hidden surface and lighting and shading problems. A large number of solutions 
exist to both of these problems. Very few solutions can be applied to both. One 
such solution is ray tracing. Ray tracing is the process of following an imaginary 
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ray from a viewpoint through a pixel on a screen and into a scene to determine if 
it intersects any objects in the scene and then calculating the intensity of the pixel 
it went through based on the the final destination of the ray. As in most cases 
where one solution is found to several problems, that solution is seldom the best 
for all the problems it is applied to and so it becomes a matter of trade-offs. Such 
is the case with ray tracing. Among the hidden surface removal techniques, ray 
tracing is the least efficient being referred to as a brute force technique. In 
contrast, it has been labeled as one of the most elegant techniques in regards to 
lighting and shading [Ref. 2: p. 137]. Because of this latter fact, ray tracing has 
become an important technique in computer graphics. Ever since the idea behind 
ray tracing was suggested by Appel, numerous articles, studies, and 
implementations have been done on it. These in turn have spawned fruther 
extensions and modifications. [Ref. 1: p. 296] 

A. DEFINITION AND OVERVIEW 

The idea behind ray tracing lies in the theory that the light in our 
environment can be modeled as rays. After being emitted from a source, the rays 
are then reflected and refracted through a scene. Some of the rays eventually find 
their way to the eye where the scene is recreated (Figure 1.1a). These light rays 
are emitted from light sources, such as the sun. An infinite number of light rays 
exist, but only a small percentage of them are received by us. To try and trace 
these rays from the source is computationally expensive. Appel suggested that 
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instead of tracing the rays from the source that they should be traced backwards 
from the viewer, thus dealing with only those rays that actually contribute to the 
scene (Figure 1.1b). [Ref. 1: p. 296] 

The basic ray tracing algorithm is a very simple one and not difficult to 
implement. The basic algorithm is a hidden surface algorithm. All hidden surface . 
algorithms can be classified based on the coordinate system or space in which they 
are implemented. These are either in object space or in image space. The ray 
tracing algorithm falls under the category of image space. This category of 
algorithm is implemented in the screen coordinate system in which the objects are 
viewed. Unfortunately, the calculations are performed only to the precision of the 
scene representation, which generally provides poor resolution. The image space 
algorithms work by comparing every object in the scene with every pixel. Such an 
algorithm is computationally expensive. Ray tracing algorithms have three parts: 
a viewpoint, a raster screen, and a set of objects (Figure 1.2). In the algorithm, 
the viewpoint is along the positive z axis. From this point, a ray is shot into the 
scene through the center of every pixel on the raster. Each of these rays is then 
traced and compared against every object in the scene to determine if there is an 
intersection with any of them. It is in the determination of a possible intersection 
point that a ray tracer spends anywhere from 75 to 95 percent of its time. If 
there is an intersection, then the intensity at the pixel is determined using the 
intersected object’s attributes and an appropriate illumination model. If there is 
no intersection, then the pixel intensity is determined by the background 
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F igure 1 . 2 

Layout of Ray Tracing Scene [Ref. 1: p.297] 
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intensity. This procedure is then repeated for every pixel on the raster. When 
the ray tracer is used as a hidden surface algorithm, intersection testing stops 
after the first intersection. Extensions to the basic ray tracing algorithm, that 
showed its usefulness in implementing a global illumination model, were originally 
implemented by Whitted [Ref. 4] and Kay [Ref. 5 and 6]. In these extensions of 
ray tracing, additional rays are calculated, specifically the reflected and refracted 
rays, and then tested to see if they intersect with any objects in the scene. This 
process of generating new rays and tracing them to check for possible intersections 
is continued until the rays either leave the scene or stack space is exceeded. In 
such a case, the remaining rays are treated as if they had left the scene. This 
process, illustrated in Figure 1.3a, for a single ray with intersections is easily 
represented using the tree structure shown in Figure 1.3b. Here each node of the 
tree represents a ray surface intersection. At each node, at least one and 
sometimes two subbranches are generated. One branch of each of the reflected 
and refracted rays is generated from the point. [Ref. 1: pp. 190-296] 

B. ORGANIZATION 

This study is broken into three areas: data requirements, ray tracing 
methodology, and the intensity problem. The first section reviews the data needed 
for a lighting and shading modeler, hereafter referred to as a Tenderer, of which a 
ray tracer is an integral part. The second section reviews the actual process of 
tracing a ray through a scene to be rendered. The third section looks briefly at the 
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Figure 1.3a 

Obtaining Global Data Reflection & Transmission Rays 
fRef. 3: p.61] 
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Figure 1.3b 

Global Data Tree [Ref. 3: p-62] 
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illumination problem and how it relates to the ray tracer, 
chapters present the implementation, and known limitations of 
with areas of future research. 



The concluding 
the model along 
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II. DATA REQUIREMENTS 



A. OVERVIEW OF THE DATA REQUIREMENTS 

The importance of the intersection routines in the ray tracer is apparent in 
the fact that a ray tracer spends 75 to 95 percent of its time determining 
intersections [Ref. 1: p. 297]. The key to determining intersections, however, lies 
in large part on the data used to describe the scene that is being rendered. 
Information is needed not only to describe the entire scene that is being rendered 
but more importantly to describe each object in the scene. Scene data is that 
information needed to completely describe a picture, i.e., the number, kind, shape, 
and color of any objects in the picture along with the background intensity and 
light source information. This information must be properly ordered and broken 
down. Falby [Ref. 3] suggested that a scene be broken into three categories: 
object, view, and light. Each of these areas is examined below in the context of a 
ray tracing algorithm. 

1. Object Data 

The data pertaining to each object in the scene can be grouped into two 
categories: polygon and bounding volume. The reason for this breakdown is 
twofold. First, each object in the scene is composed of polygons. They are the 
basic building blocks of the scene. Second, in order to reduce the number of 
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intersection checks, it is necessary to set up a boundary around each object so 
that the ray tracer only performs the intersection checks in the regions that 
actually contain an object. Such a boundary is called a bounding volume. We 
examine the polygon data first, 
a. Polygon Data 

Since the main focus of the ray tracing algorithm lies in determining 
the intersections between the rays shot into a scene and the objects that make up 
a scene, and since each object is comprised of polygons, the problem is really one 
of determining the intersection points between the rays and the polygons. From 
the fundamentals of vector calculus, it is known that in order to determine the 
intersection between a ray and a polygon only the vertices of the polygon are 
needed as well as the direction of the ray and a point on the ray. Since the object 
is to be constructed of polygons, its vertices are known. Therefore, it is only 
necessary to ensure that these points are stored in some manner, such as a record, 
so as to be accessible to the ray tracer. In order to determine what the intensity 
of the pixel is through which the ray passes, it is essential that the characteristics 
of the object whose polygon was intersected be available. Since an object is made 
of polygons, they inherit the characteristics of the object. These characteristics 
also need to be readily accessible and, therefore, need to be stored in some 
manner. The following is a list of the basic object characteristics that need to be 
available: (l) the specular, diffuse, and transmission coefficients; (2) the Phong 
specular exponent; and (3) the index of refraction, see Table 2.1. [Ref. 3: p. 68] 
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TABLE 2.1: OBJECT DATA 



FIELD NAME 


VARIABLE NAME 


VALUE 


Polygon Vertices 


x, y, 2 


real 


Diffuse Coefficient 


Kd , 

r a o 


real (0-1) 


Specular Coefficient 


ijr 


real (0-1) 


Transmission Coefficient 


,g ,6 


real (0-1) 


Unit Surface Normal 


x, y, 2 


real (0-1) 


Phong Specular Exponent 


n 


integer (0-200) 


Index of Refraction 


^2 


real 



b. Bounding Volume Data 

The major disadvantage of ray tracing is that it takes so much time. 
This is hard to avoid since it is so computationally expensive. It is essential, 
therefore, that more efficient techniques be developed to assist in reducing the 
number of calculations. Several techniques already exist with the bounding 
volume being the most effective [Ref. 1: p. 298]. In the description of ray tracing 
given so far, it has been stated that a ray is checked to see if it intersects with any 
object. Upon dissecting this statement further, a better understanding of the 
intersection problem can be realized. Unless some optimization is done, the ray 
tracing algorithm is forced to do the following. Each ray must be checked for a 
possible intersection with each object. Since each object is made up of polygons, 
then the ray must be checked for a possible intersection with each polygon. For a 
complicated object, such as a teapot, xhis requires a large number of checks and 
must be done for each object. The purpose for establishing the bounding volume 
lies in two facts. The first is that generally scenes are mostly background with 
just a few objects, hence very few of the rays actually hit anything. Therefore, 
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most of the intersection tests done are a waste of time. Second, a ray can only hit 
one object at a time. To have it process through the entire list of objects, when 
intersections with most of them can be eliminated, is needless. A bounding 
volume is, therefore, a method of enclosing each object in the scene in a simple 
containment vessel, which in effect creates a boundary around the object. Once 
this boundary is established, the number of overall intersection tests can be 
greatly reduced, as in the example of a teapot, which might easily have over a 
hundred polygons. If it is surrounded by a bounding box consisting of just six 
polygons, the number of intersection tests can be significantly reduced. In this 
situation, instead of having to test each ray against each polygon of the object, 
only those rays that penetrate the bounding volume need to be checked. Thus 
the bounding volume is a way to filter out unnecessary intersection tests by 
limiting the tests to those rays that Eire most likely to intersect an object. 

Just as the use of a bounding volume greatly increases the efficiency 
of the ray tracing algorithm, the use of the right kind of a bounding volume can 
improve upon that even more. In Rogers [Ref. l], the bounding volumes 
suggested are a bounding box and a bounding sphere (Figure 2.1), each of which 
has advantages and disadvantages. The bounding sphere is much easier to 
implement although it is less efficient in reducing the target area than is the 
bounding box, see Figure 2.1. The bounding box, on the other hand, is 
computationally expensive to implement. The data needed to establish a 
bounding sphere is minimal. It only requires a center point for the object and a 
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Figure 2.1a 

Object Surrounded by Bounding Sphere fRef . 1: p . 2981 




Figure 2 . lb 

Object Surrounded by Bounding Box [Ref. 1: p.298] 
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radius that encompasses every point of the object. The bounding box, on the 
other hand, requires far more data in that the polygons that make up the box 
must be described. 

2. View Data 

In rendering any scene, certain information can be applied to the scene 
as a whole. This information is grouped together to form the view data. This 
data consists of the viewpoint position, a constant to prevent division by zero, a 
refraction index for the global medium, the ambient light intensity, the 
background color, and the scene dimensions, see Table 2.2. [Ref. 3: pp. 74-75] 

3. Light Data 

To support a lighting and shading model, it is necessary to include 
certain information on the light source for the scene. That information must 
include the position of the light source, its intensity, its type, geometry, and 
dimension, see Table 2.3. 



TABLE 2.2: VIEW DATA 



FIELD NAME 


VARIABLE NAME 


VALUES 


Viewpoint 


x, y, 3 


real 


No Zero Constant 


Ko 


real (<0) 


Global Refraction Index 




real 


Ambient Light 


^ a r,g,b 


real (O-l) 


Background Light 


Ib r, g ,b 


real (O-l) 


Scene Size 


x, y 


integers 
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TABLE 2.3: LIGHT DATA 



FIELD NAME 


VARIABLE NAME 


VALUE 


Light Position 


y » 2 


real 


Intensity 


^r.g .b 


real (0-1) 


Type 


type 


enumerated (point distributed) 


Shape 


shape 


enumerated (circular, rectangular) 


Dimensions 


x, y 


real 



B. DATA STRUCTURE FOR A RAY TRACER 

Falby [Ref. 3] suggested a data structure for a multi-illumination model 
renderer. That data structure, with minor variations, has been used in this study. 
In [Ref. 3] a complete derivation of the data is presented, so for the purposes of 
this study only a brief description is given here. Figure 2.2 illustrates the layout 
of the data structure as used in this study. This data structure essentially consists 
of arrays of records layed out in a hierarchical organization. Starting from the 
highest level it consists of the following: a picture record, an array of light records, 
an array of objects, an array of subobjects, an array of common part records, an 
array of polygons, and three arrays for the vertices. Each of these is now 
examined. 

1. Picture 

Picture is a single record which contains the view data mentioned earlier. 

2. Lights 

The lights array is an array of records, with one record for each light 
source in the scene. Each record contains the light data mentioned above. 
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3. Objects 



The objects array is an array of records, with one record for each object 
in the scene. In this study, an object is the highest order item in a scene. Just as 
the scene is divided up into objects, each with its own bounding volume, so is 
each object broken down into subobjects, each with its own bounding volume. 

4. Subobjects 

The subobjects array is an array of records, with each array belonging to 
one object record. For example Figure 2.3a shows one object, a barbell, that is 
divided into three subobjects which are: the left weight, the right weight, and the 
bar. The record layout for this is as illustrated in Figure 2.3c. A subobject is the 
smallest item in the scene. Each object has at least one subobject. A subobject is 
composed of polygons or it is a sphere. Using Figure 2.3a as an example again, 
the left and right weights are spheres and the bar, instead of being a perfect 
cylinder, is composed of polygons and actually has an octagonal shape, Figure 
2.3b. Aside from containing a pointer to the common part record, examined next, 
and data for its bounding volume, it also contains information on the subobject 
type. This subobject type field indicates the geometry of the subobject, i.e., it is 
either a sphere or a polygonal object, which is an object composed of polygons. 
This information is stored because different, intersection routines are used for each 
object type. Currently a 1 indicates that planar intersection routines should be 
used and a 0 indicates spherical intersection routines should be used. 
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5. Common Parts 



A good illustration of common parts is found in an ordinary 
checkerboard. Figure 2.4a. In this figure one object exists— the checkerboard. This 
in turn has one subobject, itself. This subobject has two common parts: the white 
squares and the black squares. Table 2.1 listed the characteristics of an object 
and it is in the common parts record that these characteristics are stored. Each of 
these common parts records contains a pointer to an array of polygon records. It 
is through this arrangement that the polygons inherit the characteristics of the 
object. Therefore, the common parts array, also called the Cparts array, is an 
array of records, with one array per subobject, and each common parts record 
points to its own set of polygons, Figure 2.4b. 

6. Polygons 

The polygons array, too, is an array of records with one array for each 
subobject. This is the smallest physical item in the scene and the one against 
which the actual intersections are determined. 

7. Vertex Array 

The vertex array is an array of points that define the polygons that 
compose the subobject. 

This data structure as presented by Falby [Ref. 3] proved itself to be both 
flexible and easy to use. An example of a data base that used this structure and 
which was used in testing this ray tracer can be seen in Appendix B. 
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Figure 2.4a - Example of Object with Two Common Parts 
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Figure 2.4b - Record Layout of 2.4a 
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III. RAY TRACING INTERSECTION CONSIDERATIONS 



The methodology behind ray tracing is quite simple. However, it does require 
an understanding of the fundamentals of vector calculus and geometric optics. A 
discussion of the fundamentals is beyond the scope of this study. An overview of 
some of the fundamentals is in order. 

A. RAY TRACING MECHANICS 

By its very definition, ray tracing is simply the tracing, or following, of a ray 
from its source through space and determining any possible intersections that may 
occur between it and an object. The natural way to model a ray in order to do 
this tracing is by using vectors. A vector is not only a precise way to represent a 
ray but the basic operations on vectors in three space, addition, subtraction, dot 
product, and cross product provide the tools necessary to determine the 
intersections. These tools, along with other techniques found in vector calculus 
and geometric, optics provide the means to deal with the two problems 
encountered in ray tracing, i.e., the ray direction determination problem ans the 
intersection problem. 

1. The Ray Direction Problem 

Solving the ray direction problem is both the first and last step 
encountered in the ray tracing process. Determining the initial ray from the view 
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position, usually referred to as the view ray, is the simplest to solve. Every point 
in a coordinate system can be associated with a ray, and determining the direction 
of a ray between two points can be solved by using vector subtraction. The last 
step in the ray tracing process is determining what takes place when a ray 
intersects an object. This requires the application of the laws of geometric optics. 
Once a ray strikes an object, either one or two additional rays will be generated. 
These new rays are referred to as the reflected and refracted rays, Figure 3.1. The 
three basic laws of reflection and refraction are listed as [Ref. 8: pp: 32-33] 

1. The incident, reflected, and transmitted rays all reside in a plane, known 
as the plane of incidence, which is normal to the surface of the object. 

2. The angle of incidence is equal to the angle of reflection © t = 0 r . 

3. The incident and transmitted ray directions are related by Snells’ law: 

TrsinOj = n t s'mQ r 



An illustration of these laws is shown in Figure 3.1. Rogers [Ref. 1: p. 367] 
provides a method for determining the direction of the reflected and refracted 
rays. The direction of r, the reflection ray, and p, the refraction ray are given as: 



where 



r = v' +2 n 

p — kj(n t ' ) — li 



v 

/ 

V = 

vn 
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Figure 3.1 



Creation of Reflected and Refracted Rays [Ref. 3s p.55] 
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k f = (* 2 I v '\ 2 -\ v ' + ”l 2 ) 2 
^2 

k » = — 

Here k is the ratio of refractive indices, k f is the Fresnel coefficient, v is the 
normal vector in the direction of the incoming ray, v' is the unit normal vector in 
the direction of the incoming ray, n is the unit surface normal, and rj ^ and tj 2 are 
the refraction coefficients for mediums the rays pass through. This is illustrated 
in Figure 3.2. 

2. The Intersection Problem 

As stated above the intersection computation is the most time 
consuming part of the ray tracing process. It is not that the process itself is so 
difficult but because several steps need to be done for each iteration. Two types 
of intersection computations are required to be performed: determining the 

intersection between a line and a sphere and determining the intersection between 
a line and a polygon. The first type is the simplest to solve and is why the sphere 
is generally used as the bounding volume. The calculation of the intersection 

point between a line and a sphere involves solving the equation for the line and 

the sphere simultaneously. The sphere is defined by the equation 

(x - af + (y - fff + [z - 6 ) 2 = r 2 (3.1) 

where (a, /?, 6) is the center point, r is the radius, and (x, y, z) is a point on the 

sphere. The line is defined by the parametric equations 
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Figure 3.2 

Direction of Reflected and Refracted Rays [Ref. 3: p.57] 
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(3.2) 



x = at + x Q y = bt + y Q z = ct + z Q 

where (x Q , y Q , z Q ) is a known point on the line, and a, b, and c are coefficients 

from the vector at + bj 4- ck which is parallel to the line. These equations must 

first be solved for t. The solution to this provides two answers. First, it indicates 

whether or not an intersection actually takes place. Second, if there is one, it 
indicates how many intersections, either one, in the case where the line is tangent 
to the sphere, or two, in the case where it actually enters the sphere. In the case 
of two intersection points, a check must be done to determine which is closer to 
the origin of the ray. 

The intersection between a line and a polygon is more involved. This 
problem is comprised of two parts: determining the intersection point between a 
line and a plane, and determining whether or not the intersection point lies within 
the polygon. Like the line-sphere intersection problem, this one also involves the 
solving of two equations simultaneously. The first of these is the equation for a 
plane which is defined as 

Ax + By + Cz = D (3.3) 

where A, B, C, and D are constants and (x, y, z) is a point on the plane. The 
second equation used is the parametric equation that defines a line. Eq. 3.2. The 
solution of this requires first substituting the equations for the line into the 
equation of the plane. The result is an equation in t, which when solved and 
substituted back into the equations for the line, provides the intersection point 
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between the line and the plane. Once this point is calculated, it is then necessary 
to determine whether or not it lies within the polygon. A general way of solving 
this is to simply determine the relationship between the intersection point and 
each edge of the polygon. The point that lies on the inside of each edge also lies 
within the polygon. If it fails the test for any edge, then it lies outside the 
polygon. The drawback to this approach, is that it only works for convex 
polygons. In this study we assume all polygons are convex. [Ref. 9]. 



B. THE RAY DATA STRUCTURE 

Rogers [Ref. l] suggests a data structure for a ray in a ray tracer. It is that 
data structure which is used in this study. Table 3.1 lists the data used to model 
each ray. We examine each item of this structure as adapted from Rogers [Ref. 1: 
p. 373]. 



TABLE 3.1 - RAY DATA 



FIELD NAME 


VARIABLE NAME 


VALUE 


Ray Type 


type 


enumerated or coded 


Ray Origin 


x, y, z 


real 


Ray Vector 


x, y, z 


real 


Source Ray Type 


Stype 


enumerated or coded 


Intersection Flag 


flag 


boolean or coded 


Object Index 


obj idx 


integer 


Subobject Index 


subobj idx 


integer 


Common Part Index 


Cpart idx 


integer 


Polygon Index 


polygon idx 


integer 


Intersection Point 


x, y , 2 


real 


Distance 


d 


real 


Transmitted Intensity 


h 


real 


Specular Intensity 




real 
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1. Ray Type 



The ray type field identifies a ray as either a view ray, a reflected ray, or 
a refracted ray. The values put in this field are generally of an enumerated type 
and consist of reflected, refracted, view, and none. 

2. Ray Origin 

The ray origin field contains the point that identifies the position from 
which the ray originated. For instance, if it is the view ray, its point of origin is 
the view position. If it is a reflected or refracted ray, its origin is the intersection 
point that it originated from. 

3. Ray Vector 

The ray vector field contains the vector heading of the ray. 

4. Source Ray Type 

The source ray type field contains the ray type of the source ray for this 
particular ray. For instance, a view ray does not have any source ray as it is the 
starting ray for the process. Hence, none is in the type field. If the view ray 
intersects an object and both a reflected and refracted ray are generated, then the 
source ray for both of them is the view ray. Likewise if the reflected or refracted 
ray hits an object and creates further rays then it becomes the source ray for 
those rays it creates. 

5. Intersection Flag 

Originally the intersection flag is set to false and it is only set to true 
when there is an intersection between this ray and an object. 
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6. Object Index 



The object index provides an index into the array of object records 
making it possible to select any object easily. 

7. Subobject Index 

The subobject index provides an index into the array of subobject 
records and helps uniquely identify each subobject. 

8. Common Part Index 

The common part index provides the index into the array of common 
part records uniquely identifying each common part record. 

9. Polygon Index 

The polygon index provides the index into the array of polygon records 
and uniquely identifies each polygon. 

10. Intersection Point 

The intersection point field holds the position of the intersection point 
between the current ray and an object. 

11. Distance 

The distance field contains the distance between the current ray’s point 
of origin and its point of intersection. 

12. Transmitted Intensity 

The transmitted intensity field contains the red, green, and blue 
intensity values, in a range between 0 and 1, of the light that is incoming along 
any refracted ray that this ray produces. 



37 



13. Specular Intensity 



The specular intensity field contains the red, green, and blue intensity 
values, in a range between 0 and 1, of the light that is incoming along any 
specular ray that this ray produces. 

C. INTERSECTION METHODOLOGY 
1. Intersecting a Planar Polygon 

a. Generating the Initial Ray 

The generation of the initial or view ray is shown in Figure 3.3. 
This ray, p — v, is calculated using vector subtraction. The two vectors used are 
the ones associated with the points for the view position, v, and the pixel through 
which the ray passes, p. The ray associated with the view position is to be 
subtracted from the ray associated with the pixel. 

b. Intersecting the Bounding Volumes 

After each ray is generated, each object in the object list is checked, 
one at a time, to determine whether or not the ray strikes any of the bounding 
containers of the objects in the scene. In our implementation, the bounding 
container is the sphere, which is reduced to a bounding circle, C, see Figure 3.4. 
This circle's radius is the same as the radius of the bounding sphere. 5, and it lies 
on the plane, P, which contains the center point, q, of the bounding sphere. The 
inverse of the incoming ray, t?, is the surface normal of this plane. It is not 
necessary to determine where on the bounding sphere a ray hits since at this stage 
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Figure 3.3 - Determining the Ray View 



TOP VIEW 




Figure 3.4 - The Bounding Circle 
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we are only interested in just finding out whether or not it hits it. Because we are 
not concerned with where the ray strikes the bounding sphere, a bounding circle is 
used. Determining the intersection with a circle requires less work than 
determining the intersection with a sphere. 

The first step in tracing the ray begins by taking each object and 
constructing a bounding circle. Once this circle is constructed, the intersection 
point of the view ray with it is calculated. Then the distance between the 
intersection point and the center point is determined and compared to the radius 
of the bounding circle to see if it falls within the circle. If it does fall within the 
circle, this indicates that the ray intersected the bounding volume. This process 
then needs to be repeated for each object. If an object’s bounding circle is hit, 
this process must then be repeated for each subobject of that object, 
c. Intersecting the Polygon 

Once a particular subobject is identified by the bounding volume 
tests, the common parts list is processed, one record at a time. Each one of these 
common parts records contains a pointer to those polygons that make up the 
object being rendered. This list of polygon records is then processed after the 
bounding volume processing is completed. The entire list needs only to be 
processed until an intersection is found. The processing involved in this is the 
most computationally expensive part of the entire ray tracing algorithm. This 
computation consists of three steps: determining the orientation of the polygon, 
calculating the intersection point between the ray and the plane that contains the 
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polygon, and determining whether or not the intersection point lies on the 
polygon. 

(1) Establishing the Orientation. First, each polygon needs to be 
checked for a correct orientation. This is a straightforward step carried out by 
calculating the angle between the surface normal of that polygon and the inverse 
of the view ray. If it is 90 degrees or greater, it is facing in the wrong direction to 
be intersected. If it is less than 90 degrees, the next step is to determine whether 
the ray intersects the particular plane that that particular polygon lies on. The 
first step in doing this is to determine the equation for the plane in which that 
polygon lies. 

(2) Intersecting a Plane. If the correct orientation exists for a 
polygon to be intersected, the next step is to define the plane containing the 
polygon of interest. The equation for a plane was given earlier as 

Ax + By + Cz = D 

where A, B, C, and D are constants and can be calculated by the following 
equation. 

A = yi( z 2 ~ h) + M z 3 - + y 3 ( 2 i ~ h) 

B = ~ ~ — 2 i) ~ z Z^ X \ ~ 

C = Xj(y 2 - y 3 ) + z 2 (y 3 - yj + x z {y x - y 2 ) 

D = -*i(y 2 2 3 ~ */3 2 2 ) “ * 2 (*/3 2 l - yi 2 3 ) - *3 (^1 2 2 - y 2 Z \) 

where (z p y p zj, (x 2 , y 2 , z 2 ), and (x 3 , y 3 , z 3 ) are points of the vertices of the 
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polygon. Once the equation for the plane is known, it must then be solved 
simultaneously with the equation of the line representing the ray. 

(3) Location of a point with respect to a polygon. Once the 
intersection point is determined, it needs to be checked to see whether or not it 
lies within the polygon. This is the most computationally expensive part of the 
process. This process, see Figure 3.5, requires that a plane, P, called the bounding 
plane, containing two vertices, for example A and B, of an edge and another 
arbitrary point, if, not on the polygon or its plane, be constructed. This must be 
done for each edge. Once the bounding plane is constructed, the point in question 
must then be checked to see whether or not it lies to the polygon’s side of the 
plane that now contains the edge of the polygon. This is done by plugging in 
another vertex of the polygon into the equation for the plane that wa 5 just 
constructed, and then plugging in the intersection point. If the results from these 
two equations have the same relationship, i.e. if the sign is the same, then the 
intersection point lies on the polygon side of the bounding plane. This check 
must then be repeated for each edge of the polygon. In order for the intersection 
point to lie on the face itself, it must be found to lie on the polygon’s side of each 
edge. At this point, it is determined whether or not the ray strikes the object. If 
the ray does strike the object, the pointers to this polygon must be stored in the 
ray structure along with the intersection point. This information is needed later 
on when determining the direction of reflected and refracted rays as well as the 
intensity of the light that is reflected from this position. 
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Figure 3.5 

Location of a point with Respect to a Polygon [Ref . 9] 
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2. Intersection of a Sphere 



The sphere is the easiest object to work with in a ray tracer. Since the 
sphere can act as its own bounding volume, the center point and radius are 
already available in the subobject record. This eliminates the need for the polygon 
array. The center point and radius are the only information necessary to model a 
sphere for a ray tracer. Determining the intersection of a line and a sphere is 
nothing more than the simultaneous solving of their equations for the variable t . 
The solution to this gives a quadratic equation in t which can then easily be 
solved. 
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IV. THE INTENSITY PROBLEM 



One of the strong points of the ray tracing algorithm is that a global 
illumination model can easily be integrated into it. In fact, ray tracing and global 
illumintion models seem to naturally complement each other. A global 
illumination model takes into account all of the light sources in the scene in 
calculating the intensity at each point. This means taking into account the 
ambient light that exists in a scene, light that comes directly from a particular 
light source(s), and light that is reflected off an object(s). It also includes the 
coefficients necessary to model the way an object reacts with light. A great deal 
of work has been done in this area. The most notable model is the Whitted 
illumination model, and it is the one that has been implemented here [Ref. 1: pp. 
365-366]. The Whitted algorithm is based on the three models shown in Figures 
4.1, 4.2, and 4.3. These models will now be examined more closely.* 

A. LOCAL ILLUMINATION MODELS 
1. Diffuse Reflection Model 

The first of these models is a perfect diffuser. Such a model is governed 

by Lambert’s cosine law. This law states that the intensity of light reflected from 
a perfect diffuser is proportional to the cosine of the angle between the light 

*The contents of this chapter are close adaptations from Rogers [Ref. l] and Falby [Ref. 3]. 
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direction and the normal to the surface. This can be expressed mathematically as 



n 

I - 7 ; fc rf cos0 0^0^ — 

2 

where / is the reflected intensity, I { is the incident intensity from a point light 
source, k d is the diffuse reflection constant, unique to each object, and 0 is the 
angle between the light direction and the surface normal, see Figure 4.1. Since the 
diffuse reflection coefficient k d varies from material to material and is also a 
function of the wavelength of the light, it is often easier to just assume it a 
constant for simple illumination models. [Ref. 2: p. 312] 

2. Specular Reflection Model 

The second model illustrates the characteristics of specular reflection 
which is directional, unlike diffuse reflection. This means that the greatest 
intensity of the specularly reflected light can only be seen if the view angle 
coincides with the reflection angle, Figure 4.2. The further off the viewing angle is 
from the reflection angle, the dimmer the intensity becomes. Because of the 
complex physical characteristics of specularly reflected light, an empirical model 
due to Bui-Tuong Phong is usually used for simple illumination models [Ref. 10]. 
This is expressed mathematically as 

f = I[W(® t ,\) cos 1 a 

where ^(0^ A), the reflection curve, gives the ratio of the specularly reflected 
light to the incident light as a function of the incidence angle © t and the 
wavelength A. Because u>(0 t , A) is such a complex function, it is frequently 
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Diffuse Reflection 
TRef. 1: p.3121 
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View 

V 



Light. Source 



position 




Figure 4.3 

Global Illumination Model 
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replaced by an aesthetically or experimentally determined constant k which then 
yields 

I = I.k cos n ot 

3 l S 

Also, n is a power that approximates the spatial distribution of the specularly 
reflected light. Typically a value of 200 for n is used to model a very shiny 
surface and a value of 10 is used for a dull surface [Ref. 3: p. 72]. [Ref. 1: pp. 313- 
315] 

3. Combined Model 

If just point sources are assumed, as in the two models just discussed, 
any object not receiving light directly from the source appears black. In order to 
properly render a scene, it is also necessary to take into account ambient light, the 
light that is reflected off other surfaces. Including a model for ambient light into 
the intensity calculations is not feasible. Ambient light represents a distributed 
light source and as such is a very complex function. Therefore it is treated as a 
constant diffuse term and linearly combined with the other terms. Also not 
included in the above model is the effect that distance has on light. It is well 

known that the farther away an object or light source is, the dimmer it gets. The 

L 

actual formuia to produce that proper attenuation is — , where the intensity of 

o 

dr 

light decreases as the square of the distance from the source increases. However it 
has been shown that linear attenuation can actually produce more realistic results. 
With these two additions the complete model now looks like: 
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/ = hK + 

a a 



d + K 



(fc^cos© + k } cos n a) 



where I a is the incident ambient light intensity, k a is the ambient diffuse 
reflection constant and K is an arbitrary constant that can be used to adjust the 
results. [Ref. 1: p. 313] 

The complete model just presented can now be modified to better fit in 
with a ray tracer. Recalling the formula for the dot product of two vectors allows 
writing the cos© els n-L and writing cosa as R-S which gives us: 

h 

1 - i.K + + M B ' S > 1 

d + K 



So far we have only examined the case where just one light source is present. If 
there are several light sources, the effects are added linearly, and the equation now 
becomes: 



m L i 

i = KK + E — —IM" •■£,•) + *,(*,•«)"] 
d + K 



This then is the complete local illumination model. [Ref. 1: pp. 312-316] 



B. GLOBAL ILLUMINATION MODEL 

The complete local illumination model just presented forms the basis for the 
algorithm that was implemented for this study [Ref. 1: pp. 363-378] see Figure 
4.3. 



49 



1 = k a 1 a + k d S 7 /j “ ’ £ ] + ^E 7 i,( § ‘ R,j + k J s + k J t 

i i 

In the above equation k a , k d , k $ , and k { are the ambient, diffuse, specular 
reflection, and transmission coefficients, all of which have values between 0 and 1. 
I a , I s , I t , and I t are the intensities of the ambient light, the specularly reflected 

light, the transmitted light, and the light directly from a light source. These also 
hold values between 0 and 1. The remaining variables n, L ., S, and axe the 
surface normal at the intersection point, the direction of the jth light source, the 
local sight vector, and the local reflection vector from the jth light source. A 
careful comparison between this model and the complete local illumination model 
reveals that the only new terms are the 7^and I t terms. These are the terms used 
to account for the light that comes in along the reflected and refracted rays that 
originated at this point. / holds the intensity for the reflected ray and I t holds 
the intensity for the refracted ray. These two values in turn are calculated using 
the exact same model. For the last intersection point in the scene, the one whose 
reflected and refracted rays do not intersect anything, 7 and I t are set to 0. The 
k } and k f terms are coefficients included to better model the way this object reacts 
with + he light incoming along the reflected and refracted rays. 

This then is the complete global illumination model used in this study. It is 
the simplicity of this algorithm that makes it so easy to understand and 
implement. In essence, it is saying that the output intensity is nothing more than 
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a sum of all possible light sources with the coefficients determining the intensity of 



light that comes from a particular object. 
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V. RAY TRACING ALGORITHM 



Rogers suggested an algorithm for ray tracing [Ref. 1: pp. 374-377]. It is that 
algorithm that has been the basis for this study. The following is a description of 
that algorithm as it has been implemented here. 

A. TRACING THE RAYS 

To begin the ray tracing process, the first thing done is the determination of 
the direction of the view ray.* The ray data, mentioned in chapter III, is then 
initialized. After the view ray is generated and the ray data initialized, the ray, 
which is represented by this ray data, is pushed onto the stack, which is used to 
implement the ray tracing tree. The process then moves into the actual ray 
tracing loop. 

Once in the loop, the stack is checked to see if it is empty. If it is not empty, 
the stack is then popped to access the ray information. The first thing checked is 
whether or not the intersection flag is set. If it is, that indicates that that 
particular ray has already been terminated (by hitting an object), and the process 

of determining the intensities begins. 



‘Each ray, is modeled as a vector, and is converted to a unit vector immediately after its 
determination. 
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If the intersection flag is not set, a new ray has been generated, either a new 
view ray or a new reflected or refracted ray. All of these rays are grouped under 
the more general title of a shooting ray. At this stage, the ray must be sent 
through the intersection procedures to determine whether or not it intersects any 
object. The intersection routines start at the highest level of the picture and step 
through the linked list of objects, subobjects, and common parts to the actual 
polygons. 

At the top level of the intersection routines, each one of the objects is 
checked. First, they are checked to see if there is an intersection with the object’s 
bounding circle. If there is an intersection, the distance between the intersection 
point and the origin of the ray is calculated and placed in the ray data. Second, a 
check is done to insure that the objects lie in front of the ray’s origin. In front 
refers to those objects that lie in the direction of the shooting ray from its point of 
origin. Since the ray is being modeled by a vector, which only indicates direction, 
every object along the line described by the vector and the origin of the ray is 
considered, see Figure 5.1. The way to test for this is to generate a test vector 
between the origin of the ray and any intersection points of the shooting ray, 
eliminating the intersection points that lie in the opposite direction from the 
shooting ray. 

As each object is processed, the objects whose bounding volumes are 
intersected by the shooting ray are processed further to see if any of their 
subobjects axe also intersected. The same basic procedure for finding an 
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Figure 5.1 — Object Locations 
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intersection is used here as is used for finding an intersection for the objects. If the 
object and the subobject are one and the same, this check need not be done. 

For each subobject whose bounding volume is intersected, the next step is to 
check the type of geometric object, i.e., a sphere or a polygonal object. This needs 
to be done in order to determine which intersection routines are needed. In the 
case of a sphere, it becomes just a matter of solving two simultaneous equations, 
discussed earlier-one for a line, the other for the sphere. When dealing with a 
polygonal object, it becomes more complicated because each face of the polygon 
must be checked. The first thing that needs to be done is to check the orientation 
of the face in question. When the polygon does have the right orientation, the 
intersection process continues. The first step is to determine the equation of the 
plane that contains the polygon. This plane is calculated from any three vertices 
of the face in question. Once the plane has been established, the intersection point 
between it and the shooting ray is calculated. With the intersection point thus 
established, the next step is to determine whether or not the point lies on the face 
of the polygon. If the intersection point is found to lie on the face of the polygon, 
then this point is saved and placed in the ray data. 

If no intersection is found and the shooting ray is either a reflecting or 
refracting ray, then it is discarded and the intensity calculation process begins. If 
the shooting ray is a view ray, then the intensity is set to the background 
intensity which is then displayed and the next ray is generated. 
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When there is an intersection, the first thing to be checked is whether or not 
there is enough room on the stack. Since the stack holds only part of the ray tree 
at any one time, it need only be long enough to contain the longest anticipated 
branch. A particular branch of the ray tree is terminated when both the reflected 
and refracted rays at an object intersection leave the scene or when the available 
stack length is exceeded. When both rays leave the scene, their contribution to 
the illumination at the source ray is zero. When the available stack length is 
exceeded, the algorithm calculates the illumination at the source ray using only 
the ambient, diffuse, and specular reflection components at the source ray 
intersection. An extension algorithm is given in Rogers whereby the algorithm can 
be extended one additional depth in the tree without exceeding the maximum 
stack depth. However, the implementation of that was not necessary here. When 
the stack does get full, it becomes a matter of calculating the intensity at that 
intersection point and setting the appropriate value I t or I $ in the source ray. 
[Ref. 1: p. 372] 

When the stack is not full, then the distance between the source point of the 
shooting ray and the intersection point is determined and placed in the ray data. 
The ray is then placed back on the stack. Once that is done, any reflecting and/or 
refracting rays emanating from the intersection point are determined, their ray 
data initialized, and then placed on the stack as new rays with the reflecting ray 
being placed on first. It is important to keep this order of rays in mind because it 
is necessary to know the number of rays to pop when setting the intensities of the 
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source rays. With these new rays in place on the stack, the program returns to 
the beginning of the ray tracing cycle. In the absence of reflecting or refracting 
rays, then the first ray popped is the view ray. Now since this ray already has its 
intersection flag set, the intensity at the point of its intersection is calculated and 
displayed. If the ray popped is a reflective or refractive ray, then it becomes the 
new shooting ray. This cycle continues until either no more reflecting or 
refracting rays are produced or until the stack becomes full. This process is 
summarized by the pseudocode description in Figure 5.2. 

B. DETERMINING THE INTENSITY 

If a ray’s intersection flag is found set at the beginning of the ray tracing 
process, that ray is sent into the intensity processing part of the ray tracer. The 
first step in this process is to take the ray data and to determine the intensity at 
that ray’s intersection point. This intensity is then divided by the distance 
between the ray’s point of origin and its intersection point in order to properly 
attenuate it. This process is done for each of the primary colors — red, green, and 
blue. If a view ray is being considered, that means that it was the last ray on the 
stack. Therefore, it is the final intensity to be calculated and hence is the actual 
one displayed. If it is a refracted ray. then the intensity just calculated becomes 
the I t intensity in the source ray for the ray currently being examined. If it is a 
reflected ray, then the intensity just calculated becomes the I $ intensity in the 
source ray for the ray currently being examined. The stack is set up so that the 
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READ IN DATA FILE 



for Y := 1 to MAX_ROWS do 
for X := 1 to MAX_COLUMNS do 
INITIALIZE VIEW RAY 
PUSH RAY ONTO STACK 
repeat 

POP RAY FROM STACK 
if INTERSECTION FLAG SET then 
CALCULATE INTENSITY 
else 

CHECK FOR INTERSECTION 
if INTERSECTION FLAG SET then 
if STACK EXCEEDED then 
CALCULATE INTENSITY 
else 

PUSH RAY BACK ON STACK 
CALCULATE REFLECTED RAY 
CALCULATE REFRACTED RAY 
if REFLECTED RAY EXISTS then 
INITIALIZE REFLECTED RAY 
PUSH REFLECTED RAY ON STACK 
end if 

if REFRACTED RAY EXISTS then 
INITIALIZE REFRACTED RAY 
PUSH REFRACTED RAY ON STACK 
end if 
end else 
else 

if (CURRENT RAY TYPE = VIEW RAY) then 
SET INTENSITY TO BACKGROUND COLOR 
end else 

until STACK EMPTY 
DISPLAY PIXEL 

end FOR-X 
end FOR-Y 



Figure 5.2 - Pseudocode Description of the Ray Tracing Process 
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source ray is on the bottom of the ray tuple with the reflected ray above it and 
the refracted ray above that. In order to set these values, it is necessary to pop 
the stack to gain access to the source ray. Once these values are set in the source 
ray, the program returns to the beginning of the loop and pops the next ray off 
the stack discarding the last ray as it is no longer needed. [Ref. 1: p. 377] 

The intensity algorithm [Ref. 1: p. 377] although very simple in design, 
depending on the number of light sources in the picture, can also become a time 
intensive part of the ray tracing program. The entire ray data set, listed in Table 
3.1, is sent into this routine as well as the pointers to the object and light source 
list. As in the intersection routine, rays are generated here. In Rogers, they are 
referred to as shadow feelers and the same term is used in this study. These 
shadow feelers are those rays represented as the vectors from the point of 
intersection to the light source, see Figure 5.3. They are used to determine the 
intensity contributed to that point from that light source. Once these rays are 
generated, they also pass through the intersection routines in order to determine 
which objects, if any, the light rays pass through en route to the intersection 
point. The first test that must be done is to determine if any of the objects passed 
through are opaque. If any are opaque, then no light reaches the intersection 
point from that light source. That point is then considered to be lying in deep 
shadow with respect to that light source. If none of the objects intersected by the 
light ray is opaque, then the light intensity needs to be attenuated accordingly for 
each of them. This attenuation entails multiplying the intensity at each point by 
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Figure 5.3 

Shadow Feelers [Ref. 3: p. 64] 
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the transmission coefficient of the object. This process then needs to be repeated 
for each light source. A running total of the intensities is maintained to be 
included in the final calculation. The calculation takes into account the ambient 
light and the light that comes in along the reflection and refraction rays. This 
process then produces the final intensity. If the input ray type is a view ray then 
it is displayed. If it is a reflected ray, it becomes the I $ value in the source ray. If 
it is a refracted ray, it becomes the I t value in the source ray. [Ref. 1: pp. 376-377] 
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VI. IMPLEMENTATION 



The prototype was written in BORLAND’S Turbo Pascal and implemented 
on an IBM AT clone. The program is 2500 lines long and takes two hours to 
generate a scene of 200x200 pixels. The scenes generated on the AT were then 
displayed using the RGB monitor of a Silicon Graphics IRIS 2400 graphics 
workstation. 

The main focus of this study was to develop a prototype ray tracer, which by 
itself is just a hidden surface removal technique. The secondary consideration was 
to integrate an illumination model into the ray tracer. Because of this focus, 
more time was spent examining the ray tracing algorithm than any of the global 
illumination models that could have been integrated along with it. 

The top three scenes in Figure 6.1 tested the ability of the prototype to 
perform as a hidden surface remover. The program proved successful in this area. 
For these scenes, a stub was used in place of the illumination model and each 
scene was lit using only ambient light. From left to right the scenes show: An 
unobstructed view of the three objects, described in a later section: the cube and 
sphere suspended above the floor but with the cube partially blocking the sphere; 
the cube and the sphere sunk part way into the floor with the cube still in front.* 

‘These were the only scenes generated to test for hidden surface removal. The remainder of 
the tests were done trying to integrate a global illumination model. 
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The testing of the global illumination model integrated into the ray tracer 
unveiled some problems. The two bottom scenes in Figure 6.1 are representative 
of the successful results. The first problem discovered was that attenuating the 
intensity by the full distance between the origin of the view ray and its 
intersection point produced totally blackened objects. The results shown in 
Figure 6.1, the bottom row, were obtained by either dividing the intensity by two, 
see the scene on the left, or by not dividing it at all, see the scene on the bottom 
right of the object. The second problem can be clearly seen by the black line that 
runs up the center of the floor on the bottom right scene. This result along with 
thosd test that generated shadows, not shown, indicated that the intensities for 
the floor were reversed. The black line, clearly seen in color Figure 5 and just 
vaguely visible in color Figure 4, is actually specular reflection and should be 
much brighter than the rest of the floor. In those scenes where shadows were 
generated the shadows were also brighter than the rest of the floor-just the 
reverse of what it should have been. 

A. INPUT 

The test data used in this study produced the scenes shown in Figure 6.1. 
This test data was in the form of a sequential file with the data structure outlined 
in Figure 2.2. The data had in it one picture record, one light in the lights array, 
and three objects in the objects array. The first object is a cube which contains 
one subobject. This subobject has one common part. The common parts record 
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then points to the six polygon records that where used to construct the cube. 
Each one of these records has its own vertex arrays. The second object is a 
sphere. This object has one subobject, and the subobject has one common part. 
Since a sphere can be its own bounding volume and since it is not constructed of 
* polygons, then the object chain for the sphere need go no further than the 
common parts array. The third object is the floor.* This consists of one 
subobject, one common part and one polygon. All objects in the scene are opaque 
and have a highly reflective surface. 

B. OUTPUT 

The output generated by the ray tracer is in the form of a bitmap, with 
values for each of the red, green, and blue components. These values range from 
0 to 1. To display these on the RGB monitor of the IRIS, each red, green, and 
blue component is then multiplied by 255 and assigned an index in the color 
table. 



‘The floor is at a 10 degree angle to the screen to provide a better perspective. 
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VII. CONCLUSIONS 



A. AREAS OF FUTURE RESEARCH 

The ray tracer is a powerful tool in computer graphics. In its original design, 
it produced the finest rendered pictures at that time. Since then there have been 
numerous extensions, among the most widely known are the ones by Phong, Blinn 
and Newell, Kay, and Whitted each of which have further enhanced the 
performance of the ray tracer [Ref. 2: pp. 343-344]. There are two main areas for 
future research: global illumination models, and intersection algorithms. Both of 
the areas are of great interest in the graphics world. Since ray tracers and global 
illumination models can be integrated so easily, working on either problem would 
undoubtedly lead to insight into the other. This would also be very easy to do 
because ray tracing naturally lends itself to a modular design making it easy to 
establish hooks for the testing of a large number of algorithms, both illumination 
and intersection. 

B. CONCLUSIONS 

We have examined the three major areas of ray tracing: the scene data 
needed, the intersection problem, and the intensity problem. The data structure 
used was adapted from [Ref. 3] and proved to be useful. The intersection 
problem, although involved, is not complex. The algorithms used in this 
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implementation are simple. The inclusion of a simple global illumination model, 
was easy to integrate and provided fair results. The ray tracer provides an 
excellent test bed program and implemented can provide a useful tool to study 
numerous problems, not only in lighting and shading but also in intersection 
determination. 
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APPENDIX A - SOURCE LISTINGS 



DECLARATIONS 



const 

{ THIS IS THE CONSTANTS DECLARATION SECTION } 
maximum size j>f_s tack = 100; 
initialjpixel j = 0; 

type 

{ TYPES DECLARATION AREA } 

vertices_array = array [1..4] of real; 

polygon_rec = record 
num vertices : integer; 

vertice_x, vertice_y, vertice j : vertices jrray; 
surface jiormaljc, surface_normal_y, surface jiormalj : real 
end; 

polygon_ptr = "polygon jrray; 

polygon jxray = array [1..6] of polygonjec; 

common _part jec = record 

K ar, K_ag, K jib : real; { AMBIENT DIFFUSE COEFFICIENT } 

K_dr, K_dg, K Jb : real; { DIRECT DIFFUSE COEFFICIENT } 

K_sr, K_sg, K_sb : real; { SPECULAR COEFFICIENT } 

R Jr, Kjg, K Jb : real; { TRANSMISSION COEFFICIENT } 
obj jefraction Jndex : real; { OBJECTS REFRACTION COEFFICIENT } 
objjphong exp : integer; { PHONG’S SPECULAR EXPONENT } 
num_polygons : integer; 
polygons : polygon_ptr 
end; 

common jpart jptr = "common _part jUTay; 

common _part_array = array [1..3] of common_part _rec; 



sub jbject jec = record 

num common parts : integer: 
common parts : common part ptr; 

sub J>sphere jadius: real; { RADIUS OF SUBOBJECTS BOUNDING SPHERE ) 
sub bsphere x : real; { CENTER OF BOUNDING SPHERE } 
sub bsphere y : real; 

sub_bsphere z : real; 

subobj type : integer; { 0: SPHERE, 1: PLANAR-POLYGON } 
end: 



sub object ptr = "sub object array; 
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sub_object__array = array [1..3] of sub_object_rec; 

object rec = record 

num_sub_objects : integer; 
sub_objects : sub object ptr; 

obj”bsphere_radius : real; { RADIUS OF OBJECTS BOUNDING SPHERE } 
obj_bsphere_x ; real; { CENTER OF BOUNDING SPHERE } 

obj_bsphere_y : real; 

obj_bsphere_z : real; 

opcode : integer { CURRENTLY NOT USED } 

end; 



object _ptr = " object _arr ay; 

object array = array [1..4] of object rec; 

light_rec = record 

I _r, I_g, I_b : real; { INTENSITY OF THE LIGHT } 

lightjc, light _y, light : real; { LIGHT POSITION } 

dimension!., dimension2 : real { NOT USED } 
end; 



light_ptr = ~light_array; 
light_array = array [1..3] of light_rec; 



picture rec = record 

num objs : integer; 



objects 

num_lights 

lights 



: object ptr; 

: integer; 

: light ptr; 



global_refraction_jndex : real; 
real; 



no zero 
ambient r 
ambient_g 
ambient b 



: real; 

: real; 
: real; 



back ground color r : real; 
background_color_g : real; 
background color b : real; 
view^position_x : real; 
view_position_y : real; 
view position z : real; 
screen_max_x : integer; 
screen_max y : integer: 

end: 



******* 4 ^********:<^**************** 4 ^******************* 4 ^**************** 



raytype = (none, view, reflected, refracted); 
colortype = (red, green, blue); 



{ RAY DATA RECORD } 
ray_ptr = "ray jrec; 
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ray_rec = record 
rayjtype 
ray ^origin _x 
ray_origin_y 
ray_origin_z 
rayvectorx 
ray vector_y 
ray vector _z 
ray stype 



: raytype; 

: real; { ORIGIN OF RAY } 

: real; 

: real; 

: real; { DIRECTION OF RAY } 

: real; 

: real; 

: raytype; { TYPE OF SOURCE RAY } 
intersection_flag : boolean; 

obj_idx : integer; { PATH DECSRIBING OBJECT INTERSECTED } 

subobj_idx : integer; 
cpart_idx : integer; 

polygon_idx : integer; 

intersect ion _x : real; { INTERSECTION POINT } 
intersect ion _y : real; 
intersect ion _z : real; 

{ DISTANCE BETWEEN RAY’S ORIGIN AND INTERSECTION POINT } 

d : real; 

{ 

INTENSITY OF LIGHT COMING IN ALONG THE REFRACTED RAY GENERATED BY THIS 
RAY 
} 

I_tr, I t g, I_tb : real; 

{ 

INTENSITY OF LIGHT COMING IN ALONG THE REFLECTED RAY GENERATED BY THIS 
RAY 
} 



I_sr, I_sg 5 Ijsb : real; 
ray_link : rayjptr 

end; { rayjptr } 

{.PA} 

| ********************************* VAR ************************************ 
var 

outfile : text; 
sysin : text; 
sysout : text; 

{ USED IN CONVERTING ALL RAYS/VECTORS TO UNIT VECTORS } 

x, y, z : real; 

unitx, unity, unitz : real; 

list : reai; 

ray_jiumber : integer; 

{ USED AS SHORTHAND BECAUSE THE OBJECT PATHS GET LONG } 
cpart jpath : common part ptr; 
subobj jpath : sub_object ptr; 

color : colortype; 
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{ USED TO IMPLEMENT THE STACK } 
ray top, 
ray_current, 
ray_next : ray_ptr; 

intensity red, 
intensity green, 
intensity_blue : real; 

pixel_x, 

pixel_y, 

pixel z : integer; 

ray generation number : integer; 

intersection_point_x, 
intersection point y, 
intersection_point_z : real; 

old intersection point x, 
old_intersection_point_y, 
old_intersection_point_z : real; 

temp_integerJL, 
temp_integer2, 
temp integer3 : integer; 

reflected ray x, 
reflected_ray y, 
reflected ray z : real; 



refracted ray x, 
refracted _ray_y, 
refracted_ray_z : real; 



surface_normal x, 
surface normal y, 
surface normal z : real; 



{ USED TO INITIALIZE RAYS } 

initial_ray_type : raytype; 

initial _rav origin x, initial rav origin v. initial rav origin j, : real: 
initial ray vector x. initial ray vector y, initial ray vpctor z : real: 
initial ray stype : raytype: 

initial intersection dag : Dooiean; 
initial_obj_idx : integer; 

initial_subobj_idx : integer; 

initial_cpart idx : integer; 

initial_polygon_idx : integer; 
initial_intersection_x, 
initial_intersection_y, 
initial intersection z : real; 
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initial d : real; 

initial I tr, initia.1 I tg, initial I tb : real; 

initial_I_sr, initial_I_sg, initial_I_sb : real; 

{ HOLDS CURRENT RAY WHILE IN ACTUAL RAY TRACING LOOP } 

current ray type : raytype: 

current_ray_ongin_x, current_ray_ongin_y, current_ray_origm_z : real; 

current ray vector_x, current_ray vector y, current_ray vector z : real; 

current_ray_stype : raytype; 

current_intersection_flag : boolean; 

current_obj_idx : integer; 

current subobj idx : integer; 

current_cpart_idx : integer; 

current_polygon_idx : integer; 

current intersection^, 

current _intersection_y, 

current_intersectionjz : real; 

current_d : real; 

current I tr. current I tg, current I tb : real; 
current I_sr, current_I_sg, current_I_sb : real; 

reflected ray : boolean; 
refracted ray : boolean; 

source_ray_num : integer; 
source_ray_type : raytype; 

{ POINTER TO PICTURE RECORD } 

picture : picture rec; 

{ 

USED TO ESTABLISH LINKED LIST OF OBJ, SUBOBJ, LIGHTS, CPARTS, AND 
POLYGONS. 

} 

light_cntr : integer; 
light_current : light ptr; 

obj_cntr : integer; 
obj curr : object ptr; 

subobj_cntr : integer; 
subobj curr : sub object ptr; 

cpart_cntr : integer; 
cpart_curr : common part ptr; 

poly_cntr : integer; 
poly next, 

poly_curr : polygon ptr; 
vertice cntr : integer; 
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{ USED WHEN READING IN DATA } 
poly_Joopjcnt : integer; 
light_loop_cnt : integer; 
cpart_loop_cnt : integer; 
subobj loop cnt : integer; 
object _loopjcnt : integer; 
vertice_loop_cnt : integer; 

intersection_X, intersection_Y, intersection_Z : real; 

{ USED TO IDENTIFY INTERSECTED OBJECTS } 

object _idx, 
subobj idx, 
cpart_idx, 

polygon_idx : integer; 
intersection_flag : boolean; 
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INTERSECTION PROCEDURES 



{ 

**************************** * INTPRCS6 PAS ******************************* 

* These are the intersection procedures used in the ray tracer program 

*x**XxXX*X***x*XX**x:*:*****X**x*******X***x*****Xxx****:X*xxX**XX5*:**XxX5*:X:**x 

} 

{ 

************************* gpfjERE INTERSECTION *************************** 

* CALLED FORM : CHECK J'ORJUBOBJ-INTERSECTION 

* CALLS TO : NONE 

* DESC : CALCULATES THE INTERSECTION POINT BETWEEN A RAY/VECTOR AND A 

* SPHERE. 

* INPUT : The centerpoint of the sphere and it’s radius. The direction of 

* the ray and a known point on the ray — which would be it’s 

* origin. 

* OUTPUT : A flag indicating whether or not there was an intersection and 

if there was the actual intersection point itself. 

*******X*******X***X*****XXX*X**X*X***X*XX*5|C****************************** 

} 

procedure sphere_intersection (Px, Py, Pz : real; 

Vx, Vy, Vz : real; 

Cx, Cy, Cz : real; 
r : real; 

var o_intersection_flag : boolean; 

var Sphere_x, Sphere_y, Sphere_z : real ); 

var 

a, b, c, tl, t2 : real; 

Xl, Yl, Zl : real; 

X2, Y2, Z2 : real; 

distancel, distance2 : real; 
radical : real; 
diffx, diffy, diffz : real; 

begin __ 

{ INITIALIZE } 

Sphere_x := 0; 

Sphere_y := 0; 

Sphere z := 0; 
o ^intersection Jlag := false: 

{ SET UP THE COMPONENTS OF THE QUADRATIC EQUATION } 
a := ( sqr(Vx) + sqr(Vy) + sqr(Vz)); 
b := ( (2 * Px * Vx) + 

(2 * Py * Vy) + 

(2 * Pz * Vz) - 
(2 * Cx * Vx) - 
(2 * Cy * Vy) - 
(2 * Cz * Vz) ); 
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c := ( sqr(Cx) 4 sqr(Cy) 4 sqr(Cz) 4 
sqr(Px) 4 sqr(Py) 4* sqr(Pz) - 

(2 * Cx * Px) - (2 * Cy * Py) - (2 * Cz * Pz) - sqr(r)); 

radical := (sqrt ( sqr(b) - (4 * a * c))); 

{ START COMPUTATIONS ON QUADRATIC EQUATION } 
if (radical < 0) then begin 

writeln(’ WARNING - imagionary number possible in SPHERE INTERSECTION’); 
o intersection flag := false 
end 

else begin 

if (radical = 0) then begin 

{ 

IF 0 THEN JUST ONE INTERSECTION POINT(THE LINE IS TANGENT TO THE SPHERE) 

} 

tl := (-b / 2 * a); { SOLVE FOR t } 

Sphere_x := Px + (Vx * tl); { CALCULATE POINT USING t } 

3phere_y := Py — (Vy * tl); 

Sphere_z := Pz 4 (Vz * tl); 
o_intersection_flag := true 
end 

else begin 

{ 

THERE WERE TWO INTERSECTION POINTS - ONE ENTERANCE POINT AND ONE EXIT 
POINT. SOLVE FOR BOTH ts’ 

} 

tl := (-b 4 (sqrt ( sqr(b) - (4 * a * c)))) / 2 * a; 
t2 := (-b - (sqrt ( sqr(b) - (4 * a * c)))) / 2 * a; 

XI := Px + (Vx * tl); 

Yl := Py + (Vy * tl); 

Zl := Pz + (Vz * tl); 

{ 

CALCULATE DISTANCE FOR BOTH INTERSECTION POINTS FROM THE POINT OF ORIGIN 

} 

distancel := (sqrt (sqr(Xl - Px) 4 sqr(Yl - Py) 4 sqr(Zl - Pz))); 

X2 := Px 4 (Vx * t2); 

Y2 := Py 4 (Vy * t2); 

Z2 : = Pz -l (Vz * t2); 



discance2 (sqrt (sqr(X2 - Px) — sqr(Y2 - Py) — sqr(Z2 - Pz))); 

{ COMPARE DISTANCES AND SELECT THE INTERSECTION POINT THAT IS CLOSER } 
if distancel < distance2 then begin 
Sphere_x := Xl; 

Sphere_y := Yl; 

Sphere_z := Zl; 
end 

else begin 
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Sphere_x := X2; 

Sphere_y := Y2; 

Sphere_z := Z2; 
end; 

o intersection_flag := true 
end; { ELSE } 

diffx := Px - Sphere_x; 
diffy := Py - Sphere_y; 
diffz := Pz - Sphere_z; 

{ 

A CHECK TO INSURE THAT THE POINT SELECTED ISN’T THE RAYS POINT OF ORIGIN 

} 

if ((diffx <= 0.0000) and (diffy <= 0.0000) and (diffz <= 0.0000)) then 
o_intersection_flag := false; 

end;{ ELSE ) 

end; {* SPHERE INTERSECTIONS *} 



{.PA} 

{ 

********************** calculate plane equation *************************** 

* CALLED FROM : FIND _INTERSECTED_POLYGON 

* CALLS TO : NONE 

* DESC : Calculates the constants of the equation of a plane when given 

* three point on the plane. 

* INPUT : Three vertices of a planar polygon. 

* OUTPUT : The A, B, C, and D constants for the equation of a plane. 
************************************************************************** 

} 

procedure calculate_plane equation(Xl, Yl, Zl : real; 

X2, Y2, Z2 : real; 

X3, Y3, Z3 : real; 
var A,B,C,D : real ); 



begin 

A := Yl * (Z2 - Z3) + Y2 * (Z3 - Zl) + Y3 * (Zl - Z2); 

B := Zl * (X2 - X3) - Z2 ‘ (X3 - Xl) + Z3 11 (Xl - X2); 

C := Xl * (Y2 - Y3) + X2 * (Y3 - Yl) - X3 * (Yl - Y2); 

D := -Xl * ((Y2 * Z3) - (Y3 * Z2)) - (X2 * ((Y3 * Zl) - (Yl * Z3))) 
- (X3 * ((Yl * Z2) - (Y2 * Zl))); 



end; { calculate_plane_equation } 



{ 
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******************* pijsjj) intersection point *********************** 

* CALLED FROM: FIND INTERSECTED JPOLYGON, 

* CHECK_FOR_SUBOBJ_INTERSECTION 

* CHECKFORINTERSECTION 

* CALLS TO : NONE 

* DESC : Calculates the intersection point between a ray /vector and 

a plane. 

* INPUT : Ray direction and a known point on the ray, i.e. , it’s 

point of origin. The constants (A,BjC,D) of the equation 

* of a plane. 

* OUTPUT : The intersection point. 

******************************************************************** 

} 

procedure find_intersection_point(i_A, i B, i C, i D : real; 

i_ray_x, i_ray_y, i_ray_z : real; 
i_source_x, i_source_y, i_source_z : real; 
var o_intersection_point_x, 
o_intersection point y, 
o^intersection point z : real); 



var 

t : real; 
begin 

{ SET UP FOR FINDING t FROM THE EQUATION FOR A LINE AND A PLANE } 
t := (i_D - ((i_A * i_source_x) + (i_B * i_source_y) + (i_C * i_source_z))) / 

((i_A * i_ray_x) + 

(i JB * i_ray_y) + 

(i_C * i_ray_z)); 

{ 

SUBSTITUTE t BACK INTO THE EQUATION FOR A LINE TO GET THE INTERSECTION 
POINT. 

} 



o intersection point x := (i ray_x * t) + i_source_x; 
o_intersection_point_y := (i_ray_y * t) + i_source_y; 
o_intersection_point_z := (i_ray_z * t) + i_source_z; 



[****<*<**X*****<*»IC*****X******** 



end; { find_intersection_point } 

{.PA} 

{ 

tXXXXXXtXXXXXXXXXXXXXXXXXXXXXX C AL EO * 1 

* CALLED FROM: FIND_INTERSECTED_POLYGON 

* CHECK_FOR_SUBOBJ_INTERSECTION 

* CHECK JFOR INTERSECTION 

* CALLS TO : NONE 

* DESC : This calculates the constants A,B,C,D of the equation of a plane 

given only a point on the plane and the surface normal of the 

* plane. 

* INPUT : A point and the surface normal of the plane whose equation you 
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* are trying to figure out. 

* OUTPUT: The constants A,B)C,D of the equation of a plane. 
************************************************************************* 

} 

procedure caleq(bsphere_x, bsphere_y, bsphere_z, 
vector _x, vector_y, vector jz : real; 
var A, B, C, D : real); 



begin 

A := vector_x; 

B := vector_y; 

C := vector_z; 

D := (vector_x * bsphere_x) -f 
(vector_y * bsphere_y) -f 
(vector_z * bsphere_z); 

end; { caleq } 



{ 

************************** POLYGON ORIENTATION **************************** 

* CALLED FORM: FINDJNTERSECTED POLYGON 

* CALLS TO : NONE 

* DESC : This checks to see if the incoming ray will hit the front face of 

* this polygon. It does this by comparing the angle between the 

* surface normal of the polygon and the inverse of the incoming ray. 

* If the angle is greater than 90 degrees then the ray is 

* approaching the back of the polygon. 

* INPUT: The direction of the incoming ray. The surface normal of hte 

* object. 

* OUTPUT: A boolean value TRUE/FALSE depending on whether or not the 

* polygon is facing the right direction. 
************************************************************************** 

} 

procedure polygon_orientation (view_vector_x, view_vector_y, view_vector_z : real; 

surface_normal_x, surface_normal_y, surface_normal_z : real; 
var o_good orientation : boolean); 



var 

cosine theta : real; 

length view vector : real; 
length_surface_normal : real; 
dot product : real; 

begin 

{ 

TAKE THE DOT PRODUCT OF THE INVERSE OF THE VIEW VECTOR AND THE SURFACE 
NORMAL OF THE POLYGON IN QUESTION. 

} 



78 



dot_product := (- view _vector_x) * surface normal x -f 
(- view _vector_y) * surface ^normal y -f 
(- view _vector_z) * surface _normal z; 

{ CALCULATE THE MAGNITUDE OF THE VECTORS } 
length _view ^vector := sqrt(sqr(- view jvector x) + 
sqr(- viewjvector y) — 
sqr(- view_vector_z)); 

length_surface_normal := sqrt(sqr( surface_normal x) -j- 
sqr( surface_normal_y) + 
sqr( surface normal z)); 

{ CALCULATE THE COSINE OF THE ANGLE BETWEEN THE RAYS } 

cosine_theta := dot_product / (length_view _vector * length_surface normal); 

if (cosine_theta > 0) then 
o_good _orientation := true 
else 

o good orientation := false; 
end; { polygon_orientation } 



{.PA} 

{ 

************************ Fpvjj) INTERSECTED POLYGON ************************* 

* CALLED FROM : CHECK FOR SUBOBJ ^INTERSECTION 

* CALLS TO : POLYGON_ORIENTATION 

* CALEQ 

* FINDJNTERSECTIONPOINT 

* CALCULATEPLANEEQUATION 

* DESC : This determines if there is an intersection between a line/ray 

and a polygon. If there is it calculates what it is. 

* INPUT: The direction of the shooting ray, a known point on that ray - 

it’s origin, and the object path identifying the subobject to 

* examine. 

* OUTPUT: A flag indicating whether or not a polygon was hit. If one was 

* hit the path identifying which one it was and the actual 

* intersection point itself. 

************************************************************************* 



} 

procedure And jntersecced polygon! i ray _x, i_ray_y, i_ray_z real; 

i source i source y, i source z : real; 

1 obj ldx, 

i_subobj_idx : integer; 
var o polygon_intersection x, 
o_polygon_intersection_y, 
o_polygon_intersection_z : real; 
var o cpart_idx, 

o_polygon_idx : integer; 
var o intersection_flag : boolean); 
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type 

{ 

THIS IS SET UP TO HANDLE 6 SIDED POLYGONS. AS EACH SIDE OF THE POLYGON 
IS TESTED TO SEE WHETHER THE INTERSECTION POINT LIES INSIDE OR OUTSIDE OF 
IT THE CORRESPONDING ELEMENT IN THE ARRAY IS SET EITHER TRUE OR FALSE. AN 
ARRAY OF ALL TRUE MEANS THAT THE INTERSECTION POINT LIES WITHIN THE 

POLYGON. 

} 

intersection array = array[l..6] of boolean; 
var 

point outside_polygon : boolean; 
intersections : intersection_array; 

cpart_cnt, 
polygon_cnt, 
vertice cnt : integer; 

intersection found, 
good orientation : boolean; 

cpart_ path ; common_part_ptr; 

polygon _ path : polygon_ptr; 

{ USED FOR ARBITRARY POINT TO DEFINE BOUNDING PLANE } 
anchor x, anchor_y, anchor z : real; 

markerD, 

check-point JD : real; 

polygonX, polygonY, polygonZ : real; 

A, B, C, D : real; 
xl, yl, zl : real; 
x2, y2, z2 : real; 
x3, y3, z3 ; real; 

diffx, difry, diffz : real; 

loop cnt : integer; 

begin 

{ INITIALIZE } 
o opart idx := 0; 

°— polygon idx := 0; 
o_polygon_intersection_x := 0.0; 

° polygon^ intersection^ y := 0.0; 
o polygon intersection z := 0.0; 
o intersection __flag := false; 
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cpart cnt := 1; 
polygon_cnt := 1; 
vertice_cnt := 1; 
intersection_found := false; 
good orientation := false; 

cpart_path := picture. objects "[i_obj_idx].sub_objects "[i_subobj _idx]. 
common_parts; 

polygon_path := cpart _path "[cpart_cntj. polygons; 



{ 

THIS CHECKS EACH COMMON PART FOR AN INTERSECTION WITH ONE OF IT’S 
POLYGONS. 

} 

repeat 

{ THIS LOOP CHECKS EACH POLYGON OF A COMMON PART FOR INTERSECTION } 

repeat 

point outside jpolygon := false; 

polygon _orientation(i__ray_x, i_fay_y, i_ray_z, 

polygon jpath "[polygon_cntj.surface_normal_x, 
polygon path "[polygon cntj. surface normal y, 
polygon_path "[polygon_cnt].surface_normal_z, 
good orientation); 

if good_orientation then begin 

caleq(polygon_path ^polygon _cnt].vertice_x[l], 
polygon_path ~[polygon_cnt].vertice_y[l], 
polygon_path * [polygon jcnt].vertice_z[l], 
polygon_path "[polygon_cnt].surface_normal_x, 
polygon_path "(polygon_cnt].surface_normal_y, 
polygon path "[polygon cnt], surf ace_normal_z, 

A,B,C,D); 



find Jntersection_point(A, B, C, D, 
i_ray_x, 
i_ray_y, 
i_ray_z, 

i source x, i source y, i source z, 
polygonX, polygonY, polygonZ); 



{ 

CHECK TO MAKE SURE YOU ARE NOT CONSIDERING THE SOURCE POINT(ORIGIN) OF 
THE RAY. 



} 

diffx := i_source_x - polygonX; 
diffy := i_source_y - polygonY; 
diffz := i_source_z - polygonZ; 

{ 

THIS SETS A FLAG IF YOU DO CONSIDER THE SAME POINT, THE CHECKS IN THE 
OTHER PROCEDURE SHOULD PREVENT THIS BUT JUST IN CASE. 
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if ((diffx <= 0.000) and (diffy <= 0.000) and (diffz <= 0.000)) 
then begin 



{ writeln(sysout,’set trip-wire’);} 
intersections! l] := false 
end 

else begin 

{ 

THIS LOOP CHECKS EACH EDGE OF THE POLYGON TO SEE IF THE INTERSECTION 
POINT LIES INSIDE OR OUTSIDE OF IT. 

} 

repeat 

{ THIS SELECTS THE FIRST VERTEX OF A POLYGON } 

xl := polygon_path "[polygon jcnt].vertice__x[vertice_cnt]; 
yl := polygon_path "[polygon jcnt].vertice_y[vertice_cnt); 
zl := polygon path "[poly gon_cntJ. vertice z[vertice cnt]; 



if (vertice_cnt = (polvgon_path ' 
1 )) 



polygon cnt i.num_vertices - 



{ 

WHEN YOU PICK THE NEXT TO LAST VERTEX YOU CAN SELECT THE NEXT 
CONSECUTIVE VERTEX TO ESTABLISH THE EDGE THROUGH WHICH YOU WANT THE 
BOUNDING PLANE TO PASS. YOU THEN MUST PICK ONE OF THE OTHER VERTICES (AND 
IT DOES NOT MAKE ANY DIFFERENCE WHICH ONE. I PICK THE FIRST ONE.) TO BE 
USED TO PUT INTO THE EQUATION OF THE PLANE THE RESULT OF WHICH IS COMPARED 
AGAINST THE RESULT THAT COMES FROM PLUGGING THE INTERSECTION POINT INTO 
THE EQUATION OF THE BOUNDING PLANE. 

} 

then begin 

x2 := polygon path '[polygon cnt). 
vertice_x[vertice_cnt + lj; 

:= po!ygon_path “(polygon_cnt). 

vertice y[ vertice cnt -t- l); 

:= polygon_path “(polygon_cnt). 
vertice_z(vertice_cnt + l); 



y2 



z2 



x3 := polygon_path '[polygon_cntj.vertice_x[l]; 
y3 := polygon_path *[polygon_cnt].vertice_y[l); 
z3 := polygon_path *[polygon_cnt).vertice_z[l) 
end 

else begin 

{ 

IF YOU DO NOT HAVE THE NEXT TO LAST EDGE THEN .JUST SELECT THE NEXT 
CONSECUTIVE VERTICE TO ESTABLISH THE EDGE FOR THE BOUNDING PLANE AND 
THE ONE AFTER THAT TO PLUG INTO THE EQUATON OF THE PLANE. 

} 

x2 := polygon __path *[polygon_cnt]. 

vertice_x[vertice_cnt + 1]; 
y2 := polygon path "[polygon cnt]. 

vertice_y[vertice_cnt + 1]; 
z2 := polygon_path "[polygon_cnt|. 
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vertice^vertice cnt -b 1 ] ; 

x3 := polygon_path "[polygon_cnt]. 

vertice_xjvertice_cnt + 2]; 
y3 := polygon_path ~[polygon_cnt]. 

vertice _y(vertice_cnt + 2]; 
z3 poiygon_3ach ^polygon_cntj. 
vertice z[vertice cnt -+■ 2] 

end; 

if (vertice_cnt = polygon_path "(polygon cntj.num vertices ) 

{ 

IF THE VERTICE SELECTED IS THE LAST ONE THEN JUST PICK THE FIRST VERTICE 
TO ESTABLISH YOUR BOUNDING EDGE AND THE SECOND VERTICE TO PLUG INTO THE 
EQUATION OF THE PLANE. 

} 

then begin 

x2 := polygon_path "[polygon_cnt].vertice_x[l]; 
y2 := polygon path "'polygon cnt]. vertice y l]; 
z2 := polygon_path ~(polygon_cntj.vertice_z[lj; 

x3 := polygon_path "[polygon_cnt].vertice_x(2]; 
y3 := polygon_path "(polygon_cnt].vertice_yj2]; 
z3 := polygon_path "[polygon_cntj.vertice_z[2j 
end; 



{ 

ESTABLISH THE ARBITRARY POINT THROUGH WHICH THE PLANE WILL PASS. 

} 

anchor_x := ((x2 + xl) / 2) + 10; 
anchor_y := ((y2 + y 1) / 2) + 10; 
anchor_z := ((z2 + zl) / 2) + 10; 

calculate_plane_equation( xl, yl, zl, 
x2, y2, z2, 

anchor_x, anchor_y, anchor_z, 

A, B, C, D); 



{ 

THE RESULT OF PLUGGING IN THE VERTICE OF THE POLYGON INTO THE EQUATION 
OF THE PLANE. 

} 

markerD := (A * x3) -r (B v3) — (C z 3); 

{ 

THE RESULT OF PLUGGING IN THE INTERSECTION POINT INTO THE EQUATION OF 
THE PLANE. 

} 

check_point_D := (A * polygonX) + 

(B * polygonY) + 

(C * polygonZ); 



if ((markerD <= -D) and (check_point_D <= -D)) then 
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{ 

IF THE RESULTS HAVE THE SAME SIGN THEN THEY BOTH LIE ON THE SAME SIDE OF 
THE BOUNDING PLANE. HENCE THE INTERSECTION POINT LIES WITHIN THE POLYGON 
WITH RESPECT TO THAT EDGE. 

} 

intersectionslvertice cntl := true 
else 

{ 

IF THE RESULTS DON’T HAVE THE SAME SIGN THEN THEY LIE ON OPPOSITE SIDES 
OF THE BOUNDING PLANE. AT THIS POINT THE INTERSECTION POINT HAS BEEN 
PROVEN TO LIE OUTSIDE THE POLYGON. 

} 

if ((markerD >= -D) and (check _point_D >= -D)) then 
intersections[verticejcnt] := true 
else 

intersections[vertice cnt] := false; 
vertice cnt := vertice cnt + 1 



until (vertice cnt > polygon path Upolygon_cnt].num vertices); 
end 
end; 



{ 

CHECK THE POLYGON INTERSECTION ARRAY TO SEE IF THE INTERSECTION POINT 
FAILED THE INSIDE TEST FOR ANY OF THE EDGES. 



} 



for loop_cnt := 1 to polygon_path *[polygon_cnt].num_vertices do 
if not(intersections[loop_cnt]) then 
point outside polygon := true; 



if point_outside_polygon then 
intersection found := false 
else 

intersection found := true; 



polygon cnt := polygon cnt -f 1; 
until ((polygon cnt > cpart path*[cpart cntj.num polygons) or 
(intersection_found)); 



cpart cnt := cpart cnt + 1; 

until ((cpart_cnt > picture. objects *[i_obj_idx].sub_objects ~[i_subobj_idx]. 
num common parts) or (intersection found)); 

{ SET UP THE OUTPUT FOR THE PROCEDURE } 

o cpart ldx := (cpart cnt - i); 
o polygon idx := (polygon cnt - l); 
o_polygon intersection_x := polygonX; 
o polygon_intersection_y := polygonY; 
o_polygon intersection_z := polygonZ; 
o intersection flag := intersection found; 



84 



end; { find_intersected_polygon } 



{.PA} 

{ 

******************** CHECK for subobj intersection ************************ 

* CALLED FROM: CHECK_FOR_INTERSETION 

* CALLS TO : CALEQ 

* FIND_INTERSECTION_POINT 

* SPHERE_INTERSECTION_POINT 

* FIND_INTERSECTION_POINT 

* DESC : Check to find out if the shooting ray intersects this 

* subobject’s bounding volume. 

* INPUT : The object whose bounding volume has been hit. 

* The direction of the shooting ray. 

* The origin of the shooting ray. 

* OUTPUT : A flag indicating whether or not there has been an intersection. 

* The intersection point -- if there is one. 

* The path to the intersected object. 

} 

procedure check for subobj intersection (i object_idx : integer; 

ij-ayjc, i_ray_y, i_ray_z : real; 
i source x, 
i source y, 
i_source_z : real; 
var o intersection x, 
o intersection y, 
o intersection _z : real; 
var o_subobj idx, 
o_cpart_idx, 
o_polygon_idx : integer; 
var o intersection flag : boolean); 



var 

closest _object : real; 

{ VECTOR BETWEEN THE RAYS ORIGIN AND INTERSECTION POINT } 

view polygon vector_x, 
view polygon _vector_y, 
view polygon vector_z : real; 

distance from ntersection, 
distance_from_view_position : real; 

subobj path : sub_object_ptr; 
obj path : object_ptr; 

cpart_cnt, 
polygon cnt, 

subobj cnt : integer; { USED TO GO INTO RESPECTIVE ARRAYS } 
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A, B, C, D : real; 



{ CONSTANTS FOR EQUATION OF A PLANE } 



{ 

INTERSECTION POINT BETWEEN THE RAY AND PLANE THE BOUNDING CIRCLE IS 
INSCRIBED ON. 

} 

bpiane intersection 
bplane_intersection_y, 
bpiane intersection_z : real; 

{ 

INTERSECTION POINT BETWEEN THE RAY AND A POLYGON. 

} 

polygonX, polygonY, polygonZ : real; 
intersection flag : boolean; 



begin 

o_subobj_idx := 0; 
o_cpart_idx := 0; 
o_polygon_idx := 0; 
o_intersection_flag := false; 
o_intersection_x := 0.0; 
o intersection y := 0.0; 
o intersection z := 0.0; 



{ 

ESTABLISH A DEFAULT DISTANCE WITH WHICH THE ACTUAL DISTANCES WIL BE 
COMPARED. 

} 

closest_object := 10000.0; 

subobj_path := picture. objects ~[i_object_idx].sub_objects; 
obj_path := picture. objects; 

{ THIS LOOP CHECKS EACH SUBOBJECT OF AN OBJECT } 
for subobj_cnt := 1 to obj_path ~[i_object_idx].num_sub_objects do begin 



{ FIRST ESTABLISH THE PLANE ON WHICH TO DRAW THE BOUNDING CIRCLE } 
caleq(subobj path "(subobj jcnt]. sub bsphere x, 
subobj path "'subobj cncj.sub bsphere y, 
subobj path~'subobj cnci.sub bsphere 2, 
i_ray x, i_ray_y, i_ray_z, 

A,B, C, D); 

{ FIND THE INTERSECTION POINT ON THAT PLANE } 

find intersection point(A,B>C,D, 

1 J* a y_ x > i_ray_y, i_ray_z, 
i_source_x, ijsource y, i source _z, 
bpiane intersection _x, 
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bplane_intersection_y, 
bplane_intersection z); 

{ 

DETERMINE THE DISTANCE BETWEEN THE INTERSECTION POINT AND THE CENTER 
POINT OF THE BOUNDING SPHERE/CIRCLE. 

} 

distance_from_intersection : = 

sqrt(sqr(subobj_path "[subobj cnt], sub bsphere x- 
bplane_intersection_x) + 
sqr(subobj_path "[subobj _cnt].sub_bsphere_y - 
bplane intersection y) ^ 
sqr(subobj_path "[subobj_cnt].sub_bsphere_z - 
bplane_intersection_z) ); 

{ IF THE DISTANCE IS LESS THAN OR EQUAL TO THEN YOU HAVE AN INTERSECTION } 
if (distance_from_intersection < = 

subobjpath "(subobj_cnt].sub_bsphere_radius) then begin 
if (subobj_path "[subobj_cnt]. subobj type = 0) then begin 

{ 

IF SUBOBJECT TYPE IS A SPHERE THEN YOU NEED TO USE THIS INTERSECTION 
PROCEDURE. 

} 

sphere intersection (i_source_x, i_source_y, i_source_z, 
i_ray_x, i_ray_y, i_ray_z, 
subobj_path "[subobj_cnt].sub_bsphere_x, 
subobj path "[subobj cntj.sub_bsphere y, 
subobj_path "[subobj_cnt].sub_bsphere_z, 
subobj_path "[subobj_cnt].sub_bsphere_radius, 
intersection flag, 
polygonX, polygonY, polygonZ); 
cpart cnt := 1; 
polygon cnt := 0; 

end 

else begin 

{ IF IT ISN’T A SPHERE THEN USE THESE } 

find_intersected_polygon(i_ray_x, i_ray_y, i_ray_z, 
i source_x, i_source_y, i_source_z, 
i object_idx, 
subobj cnt. 

polygonX. polygonY, poiygonZ, 
cpart cnt, 
polygon_cnt, 
intersection_flag) ; 

end; { * ELSE * } 
if intersection flag then begin 

{ 

IF THERE HAS BEEN AN INTERSECTION THE ESTABLISH THE VECTOR BETWEEN THE 
ORIGIN OF THE RAY AND THE INTERSECTION POINT. 
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view_polygon_vector_x := i_source_x - polygonX; 
view polygon_vector_y := i_source_y - polygonY; 
view polygon vector z := i source z - polygonZ; 



{ 

CALCULATE THE MAGNITUDE (DISTANCE) OF THE RAY. 

} ... . 

distance_from_viewjposition := sqrt(sqr(view_polygon_vector_x) + 

sqr(view polygon vector y) + 

sqr(view_polygon_vector_z)); 

if distance from view position < closest_object then begin 

{ 

COMPARE IT AGAINST THE OTHER DISTANCES AND SELECT THE CLOSEST ONE 

} . .. ' 
closest object := distance from view position; 

ojsubobj_idx := subobj_cnt; 

o cpart idx := cpart cnt; 

o polygon idx := polygon_cnt; 

o_intersection_flag := intersection_flag; 

o_intersection_x := polygonX; 

o_intersection_y := polygonY; 

o intersection z ~ polygonZ 

end 

end 

end 

end; {* FOR *} 

end; { check for_object intersection } 



{.PA} 

{ 



************************ FOR INTERSECTION *************************** 

* CALLED FROM: MAIN AND CALCULATE INTENSITY 

* CALLS TO : CALEQ 

* FIND_INTERSECTION_POINT 

* CHECK FOR SUBOBJECT INTERSECTION 

* DESC : DETERMINES IF THERE IS AN INTERSECTION BETWEEN THE SHOOTING RAY 

* AND THE BOUNDING VOLUME OF AN OBJECT. 



* INPUT : Direction of the shooting ray. 

Origin of the shooting ray. 

Pointer into the object array. 

* OUTPUT : Flag indicating whether or not there was an intersection. 

If there was an intersection then then the intersection point. 

* The path to the intersected polygon. 
************************************************************************* 



} 



procedure check_for_intersection (i_ray_x, i_ray_y, i_ray_z : real; 

i_source x, i_source y, i source_z : real; 
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i_picture_objects : object ptr; 
var o_intersection x, 
o_intersection y, 
o_intersection_z : real; 
var o_object_idx, 
o_subobj jdx. 
o_cpart_idx, 
o_poly_idx : integer; 
var o_intersection_flag : boolean); 



var 

closest ^object : real; 

object _cnt, 
object_idx, 
subobj_cnt, 
cpart_cnt, 

polygon_cnt : integer; 

dist an ce_from _view posit ion , 
distance from intersection : real; 



{ 

INTERSECTION POINT BETWEEN SHOOTING RAY AND THE PLANE THE BOUNDING CIRCLE 
IS INSCRIBED ON. 

} 

bplane intersection x, 
bplane_intersection_y, 
bplane_intersection_z : real; 

view bplane vector x, 
view_bplane_vector_y, 
view bplane_vector_z : real; 



{ 

INTERSECTION POINT BETWEEN THE SHOOTING RAY AND A POLYGON. 

} 

polygonX, polygonY, polygonZ : real; 

{ 

CONSTANTS FOR THE EQUATION OF A PLANE. 

> 

A, B. C. D : real: 
found intersection : boolean; 
distance : real; 
begin 

o_object_idx := 0; 
o subobj idx := 0; 
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o cpart_idx 0; 
o poly idx := 0; 
o_intersection_x := 0.0; 
o^intersection_y := 0.0; 
o intersection z := 0.0; 
o intersection flag := false; 
closest object := 10000.0; 

found_intersection := false; 

object cnt := 1; 

{ 

LOOP TO CHECK EACH OBJECT IN THE OBJECT ARRAY. 

} 

repeat 

if (((picture. objects "[object cnt].obj_bsphere_z > i source_z) and 
(i_ray_z > 0 )) or 

( (picture. objects "[object_cntj.obj_bspnere_z < i_source_z) and 
(i ray z < 0 ))) then begin 



{ 

ESTABLISH PLANE ON WHICH TO DRAW BOUNDING CIRCLE. 

} 

caleq(picture.objects "[object_cnt].obj_bsphere_x, 
picture.objects " [object _cnt]. obj _bsphere_y, 
pic t ure . obj ec t s " [ obj ec t cn t ] . o bj _bsph ere_z , 
i ray x, i_ra.y_y, i_ray_z, 

A, B, C, D); 

{ 

FIND INTERSECTION POINT BETWEEN THAT PLANE AND THE SHOOTING RAY. 

} 

findintersection_point(A,B)C,D, 

i ray x, i_ray_y, i_ray_z, 
i_source_x, i_source_y, i_source_z, 
bplane intersection x, 
bplane_intersection_y, 
bplane_intersection_z) ; 



{ 

DETERMINE DISTANCE BETWEEN THE CENTER OF THE CIRCLE AND THE 
INTERSECTION POINT 

} 

distance from intersection : = 

sqrt(sqr(picture.objects "[object cnt]. obj bsphere x- 
bplane intersection^) + 
sqr(picture. objects "[object cnt]. obj bsphere y- 
bplane intersection y) + 
sqr(picture. objects "[object cnt]. obj bsphere z- 
bplane intersection z) ); 
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{ 

IF INTERSECTION POINT LIES WITHIN CIRCLE THEN START CHECKING THE 
SUBOBJECTS THAT MAKE UP THE OBJECT. 

} 

if distancejfrom_intersection <= 

picture. objects "fobject _cntl.obj_bsphere radius then begin 

check_for_subobj_intersection (object_cnt, 

i_ray_x, i_ray_y, i_ray z, 

i_source_x, i_source_y, i_source_z, 

polygonX, polygonY, polygonZ, 

subobj cnt, 

cpart cnt, 

polygon_cnt, 

found intersection); 



{ 

DETERMINE THE DISTANCE BETWEEN THE ORIGIN AND INTERSECTION POINT OF 
OF THE RAY 
} 

if found_intersection then begin 

view bplane_vector x := i source x - bplane intersection x; 

view bplane vector y := i source y - bplane intersection y; 

view_bplane vector z := i source z - bplane intersection z; 

distance_from_viewposition := sqrt(sqr(view_bplane_vector_x) + 

sqr(view_bplane_vector_y) + 
sqr(view bplane vector z) ) ; 



{ 

SELECT ONE CLOSEST TO RAY’S ORIGIN. 

} 

if distance from viewposition < closest object then begin 
closest object := distance from viewposition; 
o object idx := object_cnt; 
o_subobj_idx := subobj jcnt; 
o cpart idx := cpart _cnt; 
o_poly_idx := polygon_cnt; 
o_intersection_x := polygonX; 
o_intersection y := polygonY; 
o_intersection_z := polygonZ; 
o intersection flag := found intersection 
end 
end 
end; 

end; 

object cnt := object cnt *+■ 1; 
until (object_cnt > picture. num_objs); 

end; { check_for_intersection } 
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INTENSITY PROCEDURES 



{ 



****************************** STACK EMPTY *************************** 

* CALLED FROM : MAIN and POP 

* CALLS TO : NONE 

* DESC : CHECKS TO SEE IF STACK EMPTY 

* INPUT : POINTER TO TOP OF STACK 

* OUTPUT : BOOLEAN VALUE - TRUE/FALSE 
************************************************************************ 



} 

function stack empty (input_ray_top : ray_ptr) :boolean; 
begin 

STACK EMPTY := input_ray_top = nil; 



end; { STACK_EMPTY } 



{ 

*************************** STACK EXCEEDED *************************** 



* CALLED FROM : MAIN 

* CALLS TO : NONE 

* DESC : CHECKS TO SEE IF STACK IS FULL 

* INPUT : CURRENT SIZE OF STACK and MAXIMUM SIZE OF STACK 

* OUTPUT : BOOLEAN VALUE TRUE/FALSE 



*********************************************************************** 



} function stack_exceeded(input_tl, 

input_t2 : integer) : boolean; 



begin 

if input_tl > input_t2 then 
stack_exceeded := true 
else 

stack exceeded := false; 



end; { STACK EXCEEDED } 



{.PA} 

{ 

************************ CALCULATE REFRACTED RAY ************************** 



* CALLED FROM : MAIN 

* CALLS TO : NONE 

* DESC : This calculates the direction of a refracted ray. 
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* INPUT . Direction of source ray, surface normal of intersected object, 

* the objects index of refraction, the global index of refraction. 

* OUTPUT : A flag set to true or false, depending on weither or not a 

* refracted ray was created. If one was then it’s direction is 

* given. 

xx************************************************************************* 

} 

procedure calculate_refracted_ray (i_ray__vector_x, 

i_ray_vector_y, 
i_ray_vector_z : real; 

i_obj_surface normal x, 
i _obj ^surface _norm al _y , 
i_obj_surface_normal_z : real; 
i_obj_ridx : real; 
i_global_ridx : real; 
var o_refracted_ray_x, 
o_refracted_ray y, 
o_refracted_ray_z : real; 
var o refracted ray flag : boolean ); 

var 

testKf, 

testKf2, 

Kf : real; 

Kn : real; 

abs_dot_product : real; 
length_of_ray : real; 
vl_x, vl_y, vl_z : real; 

begin 

{ 

THIS PRODUCES THE ABSOLUTE VALUE OF THE DOT PRODUCT OF THE INCOMING RAY 
AND THE SURFACE NORMAL OF THE INTERSECTED SURFACE. 

} 

abs dot product := Abs((i obj surface normal x * i_ray vector_x) + 

(i_obj_surface_normal_y * i_ray_vector__y) + 

(i obj surface _normal_z * i ray_vector_z)); 



if abs_dot_product = 0 then begin 

{ 

THIS IS JUST A PRECAUTION. 

} 



o refracted ray_iiag .= false; 
o refracted_ray_x := 0; 
o_refracted_ray_y := 0; 
o refracted _r ay _z := 0 
end 

else begin 

{ 

THIS PRODUCES THE UNIT NORMAL VECTOR IN THE DIRECTION OF THE INCOMING 



93 



RAY. 

} 

vl_x := i ray_vector_x / abs_dot_product; 
v 1 y := i ray_vector_y / abs_dot jproduct; 
vl z := i_ray_vector_z / abs_dot jproduct; 

length_of_ray (sqrt(sqr(vi jc) -r sqr(vi_y) -h sqr(vL_z))); 

{ THIS GIVES THE RATIO OF THE REFRACTIVE INDICES } 

Kn := i_obj_ridx / i_global_ridx; 

{ 

THE CALCULATION OF THE FRESNEL COEFFICIENT IS DIVIDED UP THIS WAY 
INORDER TO CHECK LATER FOR AN IMAGINARY DENOMINATOR WHICH INDICATES 
TOTAL INTERNAL REFLECTION. 

} 

testKf := (sqr(Kn) * sqr(length_of_ray)); 
testKf2 := sqr(sqrt(sqr(i jDbj_surface_normal_x + vl jc) H- 
sqr(i_obj _surface_normal_y -f vl y) -r 
sqr(i_obj_surface_normal J + Vl_zj)); 



{ 



if (testKf - testKf2 <= 0 ) then begin 

IMAGINARY DENOMINATOR - TOTAL INTERNAL REFLECTION IS OCCURING 



o_refracted_ray_flag := false; 

o _refracted_ray_x := 0; 

o_refracted_ray_y := 0; 

o refracted ray z := 0 
end 

else begin 

o__refracted_ray_flag := true; 

Kf := 1 / sqrt(testKf - testKf2); 

o_refracted_ray_x := (Kf * (i _obj_surface_normal_x H- vl_x)) - 
i _obj _surf ace_normal_x; 

o refracted ray y := (Kf * (i obj surface normal y + vl y)) - 
i obj surface normal y; 

o_refracted_ray_z := (Kf * (i_obj_surface_normal_z + vl_z)) - 
i_obj_surface normal_z 



end 

end; 



} 



end: { calculate refracted ray } 

{.PA} 

{ 

*********************** CALCULATE REFLECTED RAY ************************** 

* CALLED FROM : MAIN and CALCULATEJNTENSITY 

* CALLS TO : NONE 

* DESC : CALCULATES THE DIRECTION OF A REFLECTED RAY 

* INPUT : Direction of the source ray and surface normal of intersected 

* object. 

* OUTPUT : Flag indicating existance of reflected ray. If one exists 
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then it’s direction is given. 

******»************»**»*************************************************** 

} 

procedure calculate_reflected_ray (ijray vector x, 

i_ray vector y, 
i ray vector z : real; 

i o bj __s urf ace _norm al _x , 
i jobj _s urf ac e _n orm al _y , 
i_obj_surface_normal_z : real; 
var o_reflected ray x, 
o_re fleeted _ray_y, 
o reflected ray z : real: 
var o_reflected_ray_flag : boolean ); 



vax 

abs_dot_product : real; 
length_of_ray : real; 
vl_x, vl_y, vl z : real; 

Rl_x, Rl_v, Rl_z : real; 

begin 

{ 

ABSOLUTE VALUE OF THE DOT PRODUCT OF THE OBJECTS SURFACE NORMAL AND 
OF THE INCOING LIGHT RAY. 

} 

abs_dot jproduct := Abs((i_obj_surface_normal_x * ijray_vector_x) + 

(i obj surface normal y * i ray vector y) + 

(i obj surface normal z * i ray vector_z)); 

vl x :=i ray vector x / abs_dot product; 
vl_y := i_ray_vector_y / abs_dot_product; 
vl_z := i_ray_vector_z / abs_dot_product; 

{ CALCULATE UNIT NORMAL VECTOR IN THE DIRECTION OF THE INCOMING RAY. } 
Rl x := vl_x + (2 * i_obj_surface_normal_x); 

Rl y :=vl y + (2 * i obj surface normal y); 

Rl_z := vl z + (2 * i obj surface normal z); 

length_of_ray := (sqrt(sqr(Rl_x) -b sqr(Rl_y) + sqr(Rl_z))); 

{ CALCULATE REFLECTED RAY } 
o reflected ray x Rl x length ^ofj-av; 
o reflected ray y Rl y length _of_rav; 

o reflected ray z := Rl_z / lengtn_of_ray; 
o_reflected_ray_flag := true; 

end; { calculate jreflectedjray } 



{.PA} 



95 



{ 

************************ \lquLATE intensity ********************** 

* CALLED FROM : 12 and 16 

* CALLS TO : CHECK_FOR_INTERSECTION and CALCULATE_REFLECTED_RAY 

* DESC : Calculates the intensity at any given point. 

* INPUT : Color component being calculated, current ray data, 

pointer to object data, and pointer to light data. 

* OUTPUT : Intensity at a given point, either to be displayed or 

* set in the appropriate source ray. 

******************************************************************* 

} 

procedure C ALCULATE ^INTENSITY ( input _color : colortype; 

input ray d : real; 

input_ambient : real; 
input_ray_I_t : real; 

input_ray_I_s : real; 

input_ray_vector_x, 
input ray vector y, 
input ray vector z : real; 
input_number_of_lights : integer; 
input_obj_K_a : real; 

input __obj_K s : real; 

input_obj_K_t : real; 

input_obj_K_d : real; 

input_obj_phong exp : integer; 
input_intersection_x : real; 
input_intersection_y : real; 
input_intersection_z : real; 
input_obj_surface_normal_x : real; 
input_obj_surface_normal_y : real; 
input_obj_surface_normal_z : real; 
input_picture object : object ptr; 
input_light_top : light _ptr; 
var io_intensity : real ); 

var 

ans : char; 

I_d : real; 

I_1 :real; 
j, i : integer; 
distance : real; 
sight ray x, 

">igJU ray y, 
sight ray z . real; 
unit_sight_x, 
unit sight y, 
unit_sight_z : real; 
light_ray_x, 
light_ray_y, 
light_ray_z : real; 
unit light 
unit light y, 
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unit light : real; 

unit reflected_x, 

unit_reflected_y, 

unit reflected z : real; 

intersection x, 

intersection jy, 

intersection z : real; 

reflected_light_ray_x, 

reflected_light_ray_y, 

reflected light ray z : real; 

source_x, 

source_v } 

source_z : real; 

check_x, 

check_y, 

check_z : real; 

obj idx, subobj_idx, cpart_idx, polygon_idx : integer; 

intersection flag : boolean; 

obj_light distance : real; 

reflected ray : boolean; 

intersected obj_Kt : real; 



begin 

{ THIS ELIMINATES THE SAME INTERSECTION POINT FROM BEING SELECTED AGAIN.} 
if input_ray_d > 0.1 then begin 
I J :=~~ 0.0; 

io intensity := 0.0; 
reflected_ray := false; 
intersection_flag := false; 

{ THIS SETS UP THE SIGHT RAY } 
sight ray x := -input _ray_vector_x; 
sight _ray_y := -input_ray_vector_y; 
sight jray_z -input_ray_vector_z; 
distance := sqrt(sqr(sight ray x) + 
sqr(sight_ray_y) + 
sqr(sight_ray_z)); 

unit_sight_x := sight_ray_x / distance; 
unit sight y := sight rav v / distance; 
unit sight z signt ray i distance; 

for i ;= 1 to (input number of lights -}- 1) do begin 

{ THIS GENERATES THE SHADOW FEELERS } 

light ray x ;= picture. lights "[ij. light x - input_intersection x; 
light ray y := picture. lights ~[i].light_y - input_intersection_y; 
light ray z := picture. lights "(i].light_z - input_intersection_z; 
distance := (sqrt(sqr(light jray_x) + 
sqr(light _ray_y) + 
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sqr(light_ray_z))); 



{ CONVERTS IT TO A UNIT VECTOR } 
unit light_x := light_ray_x / distance; 
unit_light_y := light ray _y / distance; 
unit light _z := light_ray z / distance; 

source_x := picture. lights ~[i].light_x; 
source y := picture. lights ~[i].light_y; 
source_z := picture. lights A [i].light_z; 

{ 

CHECK TO SEE IF ANY OF THE SHADOW FEELERS INTERSECT ANYTHING. 

} 

check for intersection(-unit_light_x, -unit_light_y, -unit_light jz, 
source_x, 
source_y, 
source_z, 

input picture object, 
intersections, 
intersection_y, 
intersection^, 

obj_idx, subobj_idx, cpart_idx, polygon_idx, 
intersection_flag); 



{ 

CHECK TO INSURE THAT THE SAME POINT IS NOT CONSIDERED AGAIN, WHICH CAN 
HAPPEN. 

} 

check_x := intersect ion _x - input_intersection_x; 
check y := intersect ion _y - input_intersection_y; 
check z := intersection _z - input_intersection_z; 



{ 

PULL THE PROPER CHARACTERISTICS OF THE OBJECT OUT TO DEAL WITH THE 
APPROPRIATE COMPONENT OF LIGHT THAT IS BEING CURRENTLY DEALT WITH 

} 

if color = red then 

I 1 := picture. lights "[i]. I r; 
intersected_obj_Kt := picture. 

objects ~[obj_idx]. 

sub objects "!subobj idx]. 

common parts depart idx|. 

K _tr; 

if color = green then 

I 1 := picture. lights ^[i]. I g; 
intersected_obj_Kt := picture. 

objects "(obj idx]. 

sub objects A [subobj_idx]. 

common_parts "[epart idx]. 

Kjg; 

if color = blue then 
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I_1 : - picture. lights ~[ij.I_b; 
intersec ted _obj_Kt := picture. 

objects *[obj_idx]. 
subjobjects ~[subobj_idx]. 
common parts "[cpart idx] . 

K_tb; 

obj_light_distance := (sqrt(sqr(intersection x - 

picture. lights ~[i]. light _x) + 
sqr(intersection y - 

picture. lights "[i]. lightly) + 
sqrfintersection z - 

picture. lights "(i].light__z))); 



{ 

IF THERE HAS BEEN AN INTERSETION AND THE POINT BEING CONSIDERED IS NOT 
THE ORIGIN OF THE RAY THEN CHECK TO SEE IF THE OBJECT IS OPAQUE. IF IT IS 
THEN RETURN TO BEGINNING OF LOOP AND CHECK NEXT SHADOW FEELER. IF IT IS 
NOT OPAQUE THEN CALCULATE INTENSITY AT THAT POINT AND CONTINUE FOLLOWING 
THE RAY TO SEE IF IT INTERSECTS ANYTHING ELSE. CONTINUE THIS LOOP 
EITHER UNTIL NO MORE OBJECTS ARE LEFT OR UNTIL AN OPAQUE ONE IS 
INTERSECTED. 

} 

{ 

THIS PART OF THE CODE HAS NEVER BEEN TESTED. I JUST CODED IT AS I 
THOUGHT IT SHOULD BE FROM THE ALGORITHM IN ROGERS BOOK PP. 377. 

QUITE FRANKLY I STILL DON’T FULLY UNDERSTAND WHAT IS SUPPOSE TO TAKE PLACE 
HERE. 

} 

if intersection flag then begin 
while ((intersection_flag) and 
((check_x > 0) or 
(check_y > 0) or 
(check_z > 0))) do begin 



intersection_flag := false; 



if not( intersected jobj^Kt = 0) then begin 
if (inputjcolor = red) then begin 

I_1 := picture. lights *[i].I_r * 
intersected _obj Kt: 

end; 

if (input color — green) then begin 



I I := picture. lights ^ (i] .1 g * 

intersected_obj_Kt; 

end; 

if (input jcolor = blue) then begin 

I 1 := picture. lights ^[i] .1 b * 

intersected obj Kt; 
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end; 



calculate ^reflected ray (input ray _vector_x, 
input_ray_vector y, 
input ray vector z, 
input obj surface_normal_x, 
in p u t _o bj jsurf; ac e_no r m al _y , 
input obj surface normal z, 
reflected_ray_x, 
reflected_ray_y, 
reflected ray z, 
reflected ray); 

distance := sqrt(sqr(reflected ray x) + 
sqr(reflected_ray_y) + 
sqr(reflected_ray_z)); 

unit reflected_x := reflected ray_x / distance; 
unit ^reflected ^y := reflected _ray_y / distance; 
unit_reflected_z ;= reflected_ray_z / distance; 



io intensity := io intensity + 

((I_l * input_obj_K_d) * 

((input_obj_surface_normal_x * unit_light_x) + 
(input _obj_surface_normal_y * unit_light_y) + 

(input jobj_surface_normal_z * unit light z))) + 

((I_l * input_obj_K_s) * 

((unit_sight_x * unit_reflected_x) + 

(unit_sight_y * unit_reflected_y) + 

(unit_sight_z * unit__reflected z) )) ; 

source_x ;= intersection_x; 
source_y := intersection^; 
source z := intersection z; 

check for intersection (-unit light x, 

-unit_light_y, 

-unit_light_z, 
source_x, 
source_y, 
source z, 

input picture object, 

intersection x, 

intersection y, 

intersection^, 

obj idx, subobj idx, 

cpart^idx, polygon idx, 

intersection_flag); 

check x := intersection^ - input intersection x; 
check_y ;= intersection y - input_intersection_y; 
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check_z := intersection z - input intersection z; 

end {if} 
end {while} 
end {if} 
else begin 

calculate reilected_ray(input ray vector x, 
input_ray_vector y, 
input ray vector z, 
input _obj_surface_normal_x, 
i n pu t _obj _s u rf ace _n orm al _y , 
input obj surface normal z, 
re fleeted _r ay _x, 
reflected_ray_y, 
reflect ed_ray_z, 
reflected_ray); 

distance := sqrt(sqr(reflected_ray_x) + 
sqr(reflected ray y) -r 
sqr(reflected_ray z) ) ; 

unit reflected x := reflected_ray x / distance; 
unit_reflected_y := reflected_ray_y / distance; 
unit_reflected_z := reflected _ray_z / distance; 

io_intensity := io_intensity 4- 

((I_l * input_obj_K_d) * 

((input _obj_surface_normal_x * unit_light_x) + 

(input obj _surface_normal_y * unit_light_y) + 

(input_obj_surface_normal_z * unit flight _z) ) ) + 

((I_l * input_obj_K_s) * 

((unit_sight_x * unit_reflected_x) + 

(unit sight_y * unit_reflected_y) + 

(unit sight z * unit reflected z))); 
end { ** ELSE **} " 
end; { ** FOR ** } 

{ 

THIS IS THE STUB TO JUST HAVE EVERY OBJECT ILLUMINATED BY AMBIENT LIGHT 
io intensity := input obj K a * input ambient; 

} " " 

| **************************** ***************** ********** | 

( 

THIS IS WHERE THE FINAL INTENSITY IS CALCULATED 
i 

io_intensity := ((input_obj_K_a * input_ambient) + 
io_intensity + 

((input_obj_K_s * input_ray_I_s) + 

(input_obj_K_t * input_ray I t) ) ) ; 

{/input_ray_d or / 2 or /l} 



{ 
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THIS IS JUST TO KEEP ALL VALUES WITHIN A RANGE WHERE THEY CAN BE DISPLAYED 



} 

if io_intensity > 1.00 then 
io_intensity := 1.00; 
if io intensity < 0.00 then 
io intensity := 0.00; 



| ********** * ************************ * *** * * * * * ****** * ** * * ** | 
end 

else begin 

{ 

IF THE INPUT RAY DISTANCE IS LESS THEN ONE THAN YOU ARE CONSIDERING THE 
SAME POINT AND HENCE THE INTENSITY THERE SHOULD BE 0. 

} 

io_intensity := 0.0 
end; 



end; { CALCULATE ^INTENSITY } 



{.PA} 

{ 

************************* display PIXEL 



***************************** 



* CALLED FROM : MAIN 

* CALLS TO : NONE 

* DESC : WRITES OUTPUT TO FILE 

* INPUT : THE VALUES FOR THE RED, GREEN, AND BLUE COMPONENTS OF LIGHT. 

* OUTPUT : NONE 



********************************************************************* 



} 

procedure DISPLAY_PIXEL ( input_intensity_red, 
input intensity green, 
input intensity blue : real; 
input_pixel__x, 
input pixel y, 
input_pixel_z : real ); 



begin 

write (outfile, input intensity red:3:2); 
write (outfile/ ’/nput^ntensityj^reeniS^); 
writeln (outfile/ ’,input_intensity_blue:3:2); 

end; { DISPLAY PIXEL } 

{.PA} 

{ 

*************************** pQp *************************************** 

* CALLED FROM : MAIN, 12, 16 

* CALLS TO : NONE 

* DESC : Removes a ray from the top of the stack. 

* INPUT : Pointer to the current top top of stack. 

* OUTPUT : The ray just popped from the stack and a pointer to the new top 



102 



top of stack. 



****X**X**X*****X*****X************************************************** 



} 

procedure POP (var output ray type 
var output_ray_origin_x 
var output ray origin y 
var output ray origin z 
var output ray vector x 
var output ray vector y 
var output_ray_vector_z 
var output ray stype 
var output intersection Sag 
var output _obj_idx 
var output subobj idx 
var output_cpart_idx 
var output_polygon_idx 
var output_intersection x 
var outputjntersection y 
var output intersection z 
var output d 



: raytype; 
: real; 

: real: 

: real; 

: real; 

: real; 

: real; 
raytype; 

: boolean; 
integer; 

: integer; 
integer; 

: integer; 

: real; 

: real; 

: real: 

: real; 



var output_I_tr, output_I_tg, output_J tb : real; 
var output_I_sr, output _I_sg, output_I_sb : real; 
var io_top : ray ptr ); 



begin 



if (stack_empty (io_top)) then begin 
writeln(’STACK UNDERFLOW ERROR’) 
end { if } 
else begin 

output _ray_type 
output j*ay_origin_x 
output j-ay_origin_y 
output_ray originjz 
output_ray_vector_x 
output _ray_vector_y 
output ray vector z 
output ray stype 



:=■ 10 top .ray_type; 

io_top ".ray_origin_x; 

ray_origin_y; 
ray origin z; 
.ray_vector_x; 
= io_top ".ray_vector_y; 
= io_top".ray vector z; 
:= io top ".ray stype; 



— io_top * 
= io_top " 
= io top' 



output intersection flag := io top ".intersection flag; 



output obj idx 
output_subobj_idx 
output cpan idx 
output polygon idx 
output_intersection x 
output intersection y 
output intersection_z 
output d 
output_I_tr 
output I tg 
output I tb 
output_I_sr 
output_I_sg 



:= io top ".obj_idx; 

:= iotop ".subobj_idx; 
io top"cpart idx: 

- io_top ".polygon idx; 
io_top " intersection x: 
:= io_top ".intersection_y; 
:= io_top ".intersection^; 
:= io top ".d; 

:= io_top ^ .1 tr; 

:= io_top " .1 g; 

:= io_top "" . I tb; 

io_top ".I_sr; 

:= io top ".I sg; 
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output I sb := io_top ".I_sb; 

ray_next := io_top ".ray_link; 

{ remove old pointer } 
dispose(io_top); 

{ set top of stack pointer to new top of stack } 
io top := ray next; 

end; { else } 

end; { POP } 



{.PA} 

{ 



PI T SH ******************************** 



* CALLED FROM: MAIN, 12, 16 

* CALLS TO : NONE 

* DESC : Places a ray on the top of the stack. 

* INPUT: The current top of the stack, and the data for a new ray 

* OUTPUT: The pointer to th new top of stack. 

*********x*»**»***************************************************** 



} 

procedure PUSH ( input_ray_type : raytype; 

input_ray_origin_x : real; 

input ray origin y : real; 

input ray origin z : real; 

input ray vector_x : real; 

input ray_vector_y : real; 

input ray_vector_z : real; 

input_ray_stype : raytype; 

input_intersection_flag : boolean; 
input_obj_idx : integer; 

input subobj idx : integer; 
input_cpart_idx : integer; 

input_polygon_idx : integer; 

input intersection x : real; 

input_intersection_y : real; 

input_intersection_z : real; 

input d : real; 

input I tr, input I_tg, input I tb : real; 

input I sr, input I sg, input I sb : real; 



var io top 



: ray ptr 



begin 



new(rav current); 
rav current", ray type 

ray_current .ray_ongin_x 
ray current ".ray origin y 
ray_current ".ray origin_z 
ray_current ".ray vector x 
ray current ".ray vector y 
ray current ".ray vector_z 
ray current ".ray jstype : 

ray current ".intersection flag 



■= input rav type; 

:= input ray origin x; 

:= input_ray origin y; 

:= input ray origin z; 

:= input_ray__vector_x; 

:= input ray vector y; 

:= input ray vector z; 

:= input_ray_stype; 

:= input intersection flag; 
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ray_current 
ray current 
ray ^current 
ray_current 
ray_current 
rav current 
ray current 
ray_current 
ray current 
ray current 
ray_current 
ray current 
ray_current 
rayjcurrent 
ray current 
io_top 



.obj idx 


:= input obj idx; 


.subobj idx 


:= input subobj idx; 


.cpart idx 


:= input cpart idx; 


.polygon idx 


:= input polygon idx; 


.intersection x 


:= input intersection x 


.intersection y 


:= input intersection v 


.intersection z 


:= input intersection z; 


.d = 


= input d; 


.1 tr 


= input I tr; 


•I tg 


= input I tg; 


.I_tb 


= input I tb; 


.1 sr 


= input I sr; 


.I_sg 


= input I sg; 


.1 sb 


= input I sb; 


.ray link 


:= io top; 



:= ray_current; 



end; { PUSH } 



MAIN 



PROGRAM RAYTRACER; 

{ 

*********************************************************************** 

** PROG : RAY. PAS 
** AUTHOR : Paul G. Smith 
** DATE : 11 May 1987 

** DESC : A ray tracing prototype with a global illumination model 
** integrated into it. 

** INPUT : A sequential scene file under PICTURE.PAS. 

** OUTPUT : A bitmap file containing the red, green, and blue color 
** color components. Their values range from 0-1 and need to 

** be converted for display on an RGB color monitor. 

*********************************************************************** 

*********************************************************************** 

} 

{ * INCLUDE FILES * } 

{$1 declare6.pas} - { DECLARATION SECTION } 

{.PA} 

{SI intprcs6.pas} - { INTERSECTION PROCEDURES } 

{.PA} 

{SI procs6.pas} _ { INTENSITY and UTILITY PROCEDURES } 

{.PA} 

{ 

*********************************** j2 ******************************** 

* CALLED FROM: MAIN 

* CALLS TO: PUSH, POP, CALCULATE ^INTENSITY 

* DESC : Calculates the light intensity at a given intersection point 

* INPUT : A complete ray data record, the pointer to the light array 

* and the pointer to the object array 

* OUTPUT : The intensity at a given intersection point. If the input 

* was a view ray then this intensity will be the intensity 

* displayed. If the input ray is a reflected ray then this 

* intensity is assigned to the I s field in the source ray. 

If the input ray is a refracted ray then this intensity is 
assigned to the I t field in the source ray. 

x:»::«:X:«:x3«:x:*:xx*x*xxac:«:»::*:*X*xx*x*:*:XxxxxXX:«*xxxae*xxx:*:*:«X***x*****XX*ac**x**:*:** 



procedure 12 (i d : real; 

i_I_tr, i_I_tg, i_I_tb : real; 

i_I_sr, i_I_sg, i_I_sb : real; 

i number oMight sources : integer; 
i ambient r, i ambient g, i ambient b : real; 
i_K_ar, i_K_ag, i_K_ab : real; 
i K sr, i Kjsg, i K sb : real; 
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i_K_tr, i_K_tg, i_K_tb : real; 
i_K_dr, i_K_dg, i_K_db : real; 
i ray_type : ray type; 
i ray top : ray ptr; 
i_light_top : light_ptr; 
i obj ptr : object ptr: 
i_obj phong exp integer; 
i_ray vector_x, 
i ray vector y, 
i_ray_ve c t° r _z : real; 

i intersection x, 
i_intersection_y. 
i_intersection_z : real; 

i surface normal_x, 
i _s urf ac e _n orm al _y , 
i_surface_normal z : real; 
var io ray_generation_number : integer; 
var o intensity_red, 
o intensity green. 
o_intensity_blue ; real ); 



var 

{ SET UP TEMPORARY AREA TO HOLD RAYS POPPED FROM STACK } 
templ_ray_type : raytype; 
tempi ray origin x, 
tempi ray origin_y, 
templ_ray origin z : real; 
templ_ray_vector_x, 
tempi jray_vector_y, 
tempi ray_vector z : real; 

tempi ray stype : raytype; 

templ_intersection_flag : boolean; 
tempi obj idx, 
tempi jsubobj_idx, 
tempi cpart idx, 
tempi polygon idx : integer; 
templ_intersection_x, 
tempi intersection y, 
tempi intersection z : real; 
templ_d : real; 

tempi I tr, tempi I tg, templ_I_tb : real; 

tempi I sr, tempi I sg, tempi I_sb : real; 

temp‘2 ray type : raytype; 

temp2_ray_ongin_x, 
temp2 jray_origin_y, 
temp2 ray origin z : real; 
temp2_ray_vector_x, 
temp2 jray_vector_y, 
temp2_ray vector z : real; 
temp2 _ray_stype : raytype; 

temp2 intersection flag ; boolean; 
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temp2_obj_idx, 

temp2 subobj_idx, 

temp2 _cpart _idx. 

temp2_polygon_idx : integer; 

temp2_intersection_x, 

temp2 intersection y, 

temp2 intersection z : real; 

temp2_d : real; 

temp2 I ^tr, temp2_I_tg, temp2 I tb : real; 

temp2_I_sr, temp2_I_sg, temp2_I_sb : real; 

begin 



{ CALCULATE INTENSITY OF THE RED COMPONENT OF LIGHT } 

calculate_intensity( red, 

i_d, 

i_ambient_r, 
i_I_tr, 
i I sr, 

i_ray_vectorjc, i^ray^vector_y, i_ray_vector_z, 

i_number_of_light_sources, 

i_ K _ar, 

i_K_sr, 

i_K_tr, 

i_K_dr, 

i obj phong exp, 
i_intersection_x, 
i_intersection_y, 
i_intersection z, 
ijsurface normal x, 
i jsurface_normal_y, 
i surface normal z, 
i_obj _ptr, 
i_light_top, 
o__in tensity _red ); 

{ CALCULATE INTENSITY OF THE GREEN COMPONENT OF LIGHT } 
calculate_intensity( green, 
i_d, 

i ambient_g, 

i_I_tg, 
i I sg, 

i ray vector x. 1 ray vector y, i rav vector z, 
i number of light sources, 

i_K_ag, 

i_ K _ S g, 

i_K_tgj 

|_K_dg, 

i _o bj _ph o n g _ex p , 
i intersection^, 
i_intersection_y, 
i intersection z, 



108 



i_surface_normal_x, 
i surface normal y, 
i_surface_normal_z, 
i_obj_ptr, 
i light_top, 

0 intensity green ); 

{ CALCULATE INTENSITY OF THE BLUE COMPONENT OF LIGHT } 
calculate intensity( blue, 
i_d, 

i_ambient_b, 

i_I_tb, 

i_I_sb, 

i_ray_vector_x, i_ray_vector_y, i_ray_vector_z, 

i_number_of_light_sources, 

i_K_ab, 

i_K_sb, 

i_K_tb, 

i_K_db, 

1 obj phong exp, 
i_intersection_x, 

i intersection_y, 
i intersection z, 
i_surface_normal_x, 
i_surface_normal_y, 
i surface_normal_z, 
i_obj_ptr, 

i flight top, 

o intensity blue ); 

if (i_ray_type = view) then begin 

{ 

THE VIEW RAY IS ALWAYS THE LAST RAY ON THE STACK WHEN IT IS POPPED THE 
INTENSITY DETERMINED FOR IT IS PASSED BACK INTO MAIN FOR DISPLAY. 

} 



{ nothing } 
end 

else begin 

if (i_ray_type = reflected) then begin 

{ 

SINCE THIS IS THE REFLECTED RAY THEN JUST ONE RAY NEEDS TO BE POPPED TO 
GAIN ACCESS TO THE SOURCE RAY. 

} 

pop (templ_ray_type, 

tempi ray origin x, templ_ray_origin_y, tempi ray origin z, 

tempi ray vector x, tempi ray vector y, tempi ray vector z, 

templ_ray_stype, 

tempi intersection flag, 

tempi obj_idx, 

templ_subobj_idx, 
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t e m p 1 _c pax t _idx , 
templ_polygon_idx, 

templ_intersection_x, tempi intersection _y, tempi intersection _z, 
templed, 

tempi I_tr, templ_I_tg, templ_I_tb, 
tempi I sr, tempi I sg, templ l sb, 
ray_top ); 

{ SET INTENSITY IN SOURCE RAY } 

templ_I_sr := o_intensity_red; 
templ_I_sg := o_intensity_green; 
templ_I_sb := o_in tensity _blue: 

{ RESTORE STACK } 

push( tempi ray type, 

templ_ray_origin_x, templ_ray_origin_y, templ_ray_origin z, 

tempi ray vector x, tempi ray vector y, tempi ray vector z, 

tempi _ray_stype, 

tempi intersection flag, 

templ_obj_idx, 

tempi subobj idx, 

tempi jcpartidx, 

tempi _polygon_idx, 

tempi intersection x, tempi intersection y, tempi intersection z, 
tempi jd, 

templ_I_tr, templ_I_tg, templ_I_tb, 
tempi I sr, tempi I sg, tempi I sb, 
ray_top ); 



end 

else begin 

{ 

SINCE THIS IS THE REFRACTED RAY TWO RAYS MUST BE POPPED TO GAIN ACCESS 
TO THE SOURCE RAY. 

} 

pop (tempi ray type, 

templ_ray_origin_x, templ_ray_origin_y, templ_ray_origin_z, 

tern pl_ray_vec tor x, templ_ray_vector_y, tempi ray_vector z, 

tempi ray stype, 

templintersection_flag, 

templ_obj_idx, 

tempi subobj idx, 

tempi opart idx, 

tempi _poiygon jdx, 

tempi intersection x, tempi intersection y, tempi intersection z, 
templ_d, 

templ_I_tr, tempi I_tg, templ_I_tb, 

templ_I sr, tempi I sg , tempi I _sb, 

ray_top ); 

pop (temp2 ray type, 

temp2 ray origin x, temp2 ray origin y, temp2 ray origin z, 
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temp2 ray_vector x, temp2 rav_vector y, temp2 ray vector_z, 

temp2_ray_stvpe, 

temp2 intersection^flag, 

temp2 obj idx, 

temp2_subobj_idx, 

temp2 cpart_idx. 

temp2 polygon idx, 

temp2 intersection x, temp2 intersection y, temp2 intersection z, 
temp2_d, 

temp2_I_tr, temp2_I_tg, temp2_I_tb, 
temp2 I sr, temp2_I_sg, temp2 I_sb, 
rav top ); 

{ SET INTENSITY IN THE SOURCE RAY. } 

temp2_I_tr := o_intensity_red; 
temp2_I_tg := o_intensity_green; 
temp2 I tb := o intensity blue; 

{ RESTORE STACK. } 

push( temp2_ray_type, 

temp2_ray_origin_x, temp2_ray_origin_y, temp2_ray_origin jz, 

temp2_ray vector x, temp2 ray_vector_y, temp2 ray vector z, 

temp2 _ray stype. 

temp2 intersection_fiag, 

temp2 obj idx, 

t e m p2 _su bo bj _idx , 

temp2_cpart_idx, 

temp2_polygon__idx, 

temp2_intersection_x, temp2_intersection_y, temp2_intersection_z, 
temp2_d, 

temp2_I_tr, temp2_I_tg, temp2_I_tb, 
temp2_I_sr, temp2_I_sg, temp2_I_sb, 
ray_top ); 

push( templ_ray_type, 

tempi ray origin x, tempi ray_origin y, tempi ray origin z, 

tempi ray_vector_x, templ_ray_vector_y, tempi ray_vector_z, 

tempi _ray _sty pe, 

tempi intersection_flag, 

templ_obj_idx, 

tempi _subobj_idx, 

tempi cpart idx. 

tempi polygon idx, 

tempi intersection x, tempi intersection y, tempi intersection z, 

tempi_d, 

templ_I_tr, templ_I_tg, templ_I_tb, 
tempi I sr, tempi I sg, templ_I_sb, 
ray_top ); 



end; 

{ 

SINCE THE RAY SENT INTO THIS PROCEDURE IS NO LONGER NEEDED IT IS 
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DISCARDED, HENCE THE RAY COUNT NEEDS TO BE DECREMENTED. 

} 

io ray generation number := io_ray_generation number - 1; 
end; 

end; 

{.PA} 

{ 

********************************* J0 ********************************** 

* CALLED FROM : MAIN 

* CALLS TO : PUSH, POP, CALCULATE ^INTENSITY 

* DESC : Used to calculate intensity at node of stack storage is 

* exceeded. 

* INPUT : Incomplete ray data record, this ray can not be continued 

because there is no room on the stack for it. Also pointers 

* to the light array and object array. 

* OUTPUT : Intensity, I_t or I_s, is set in source ray. 

*******W*******Y****Y******V****************************************** 

} 

procedure 16 (i d : real; 

i_I_tr, i_Ijtg, i_I_tb : real; 

i I_sr, i Ijsg, i I sb : real; 
i number of light sources : integer; 
i_ambient jr, i_ambient_g, i_ambient_b : real; 
i_K_ar, i_K_ag, i_K_ab : real; 
i_K_sr, i_K_sg, i_K_sb : real; 
i K tr, i K tg, i K tb : real; 
i_K_dr, i_K_dg, i_K_db : real; 
i_ray_type : raytype; 

i ray top : ray ptr; 

i_light_top : light_ptr; 
i obj ptr : object_ptr; 
i_obj_phong_exp : integer, 
i_ray_vector_x, 
i_ray_vect or _y, 
i ray vector z : real; 
i intersection x, 
i_intersection_y, 
i_intersection_z : real; 

i surface normal x, 

i surface normal y, 

i surface normal z : real: 

var o intensity red, 
o intensity green, 
o intensity blue : real ); 



var 

{ STORAGE FOR POPPED RAYS } 

tempi ray type : raytype; 
tempi ray_origin_x, 
tempi _ray_origin_y, 
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templ_ray_origin_z : real; 
tempi ray vector x, 
templ_ray_vector_y, 
templ_ray_vector_z : real; 
tempi ray stype : raytype; 
tempi intersection _flag : boolean: 
tempi obj idx, 
templ_subobj_idx, 
tempi cpart_idx 5 
tempi polygon idx : integer; 
templ_intersection_x, 
temp l_intersection_y , 
templ_intersection_z : real; 
templ_d : real; 

templ_I_tr, templ_I_tg, templ_I_tb : real; 
templ_I_sr, templ_I_sg, templ_I_sb : real; 

temp2 _ray_type : raytype; 
temp*2 ray origin x, 
temp2 ray_origin_y, 
temp2_ray_origin_z : real; 
temp2_ray_vector_x, 
temp2_ray vector y, 
temp2 _ray_vector_z : real; 
temp2_ray stype : raytype; 
temp2 intersection_flag : boolean; 
temp2_obj_idx, 
temp2_subobj_idx, 
temp2 cpart idx, 
temp2_polygon_idx : integer; 
temp2 intersection_x, 
temp2 intersection y, 
temp2_intersection_z : real; 
temp2_d : real; 

temp2_I_tr, temp2_I_tg, temp2_I_tb : real; 
temp2_I_sr, temp2_I_sg, temp2_I_sb : real; 

begin 

{*** COULD INSERT TREE EXTENSION PROCEDURE ***} 

f CALCULATE RED COMPONENT OF LIGHT } 

calculate intensitv( red. 

id, 

1 ambient r, 

i_I_tr, 

i_I_sr, 

i ray vector x, i ray vector y, i ray vector z, 

i_number_of_light_sources, 

i_K_ar, 

i K sr, 

i K tr, 
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i_K_dr, 

i jo bj _ph o n g __exp , 

i_intersection_x, 

i_intersection_y, 

i intersection z, 

i_surface_normal_x, 

i_surface_normal_v, 

i_surface normal z, 

i_obj_ptr, 

i_light_top, 

0 intensity red ); 

{ CALCULATE GREEN COMPONENT OF LIGHT } 
calculate intensity( green, 
i_d, 

i_ambient_g, 

i_I_tg, 

i_I_sg, 

i_ray_vector_x, i_ray_vector_y, iray 

i_number_of_light_sources, 

i_K_ag, 

i_ K _ S g, 

i_K_tg, 

L K _ dg > 

i_obj_phong_exp, 

i_intersection_x, 

i_intersection_y, 

1 intersection z, 
i_surface normal x, 
i_surface_normal_y, 
i_surface normal z, 
i_obj_ptr, 
i_light_top, 
o_intensity_green) ; 

{ CALCULATE BLUE COMPONENT OF LIGHT } 
calculate_intensity( blue, 
i_d, 

i ambient b, 

ij_tb, 

i_I_sb, 

i ray vector x, i ray vector y, i ray 
i number of light sources, 
i“Kjib. 

ij\_sb, 

i_K_tb, 

iK_db, 

i_obj_phong exp, 
i_intersection_x, 
i_intersection y, 

1 intersection z, 
i_surface normal x, 



vector z. 



vector z, 
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i surface _normal_y, 
i surface normal z, 
i_obj_ptr, 
i_light_top, 
o intensity blue); 



if (i_rav_type = reflected) then begin 

{ 

IF INPUT RAY IS A REFLECTED RAY THEN ONE RAY MUST BE POPPED TO GAIN 
ACCESS TO IT’S SOURCE RAY. 

} 

pop (tempi ray type, 

tempi _r ay _ori gin _x, cempl_ray_origin_y, templ_ray_origin_z, 

templ_ray_vector_x, templ_ray_vector_y, tempi jray vector z, 

tempi _ray_stype, 

templ_intersection_flag, 

tempi obj idx, 

tempi subobj idx. 

templ_cpart_idx, 

templ_polygon_idx, 

tempi _intersection_x, tempi _intersection_y, templ_intersection_z, 
tempi d, 

templ_I_tr, templ_i_ig, templ_I_tb, 
templ_I_sr, templ l sg, templlsb, 
i_ray_t°p ); 

{ SET I s IN SOURCE RAY } 
tempi I sr := o intensity red; 
templ_I_sg := o_intensity _green; 
templ_I_sb := o_intensity_blue; 

{ RETORE STACK. } 

push( tempi ray type, 

tempi ray origin x, tempi ray origin y, tempi ray origin z, 

templ_ray_vector_x, tempi _ray_vector_y, templ_ray_vector_z, 

tempi _ray_s type, 

templ_intersection_flag, 

templ_obj_idx, 

tempi subobj idx, 

tempi _cpart_idx, 

tempi polygon Jdx. 

temol intersection x, tempi intersection y, tempi intersection z, 
tempi ji, 

tempi i tr, templ_i_tg, templ_i_tb, 
tempi I sr, tempi I_sg, templ_I_sb, 
ray_top ); 



end 

else begin 

{ 
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IF INPUT RAY IS A REFRACTED RAY THEN TWO RAYS MUST BE POPPED FROM THE 
STACK TO GAIN ACCESS TO IT'S SOURCE RAY. 

} 

pop (templ_ray_type, 

tempi ray origin x, tempi ray origin y, tempi ray origin z, 

tempi ray vector _x. tempi ray vector y, tempi ray _vector_z, 

temp l_ray jstype, 

tempi intersection flag, 

templ_obj_idx, 

templ_subobj_idx, 

tempi cpart_idx, 

templ_polygon_idx, 

templ_intersection_x, templ_intersection_y, tempi intersection _z, 
templ_d, 

templ_I_tr, templi_tg, templ_I_tb, 
templ_I_sr, templ_I_sg, templ l sb, 
ray^top ); 

pop (temp2 ray type. 

temp‘2 ray_origin_x, temp2_ray_origin_y, temp2_ray_origin_z, 

temp2 ray vector x, temp2 ray vector y, temp2 ray vector z, 

temp2 ray_stype, 

temp2 intersection flag, 

temp2_obj_idx, 

temp2 subobj idx, 

temp2 cpart idx, 

temp2_polygon__idx, 

temp2 intersection x, temp2 intersection y, temp2 intersection z, 
temp2_d, 

temp2_Iir, temp2_I_tg, temp2_I_tb, 
temp2 I_sr, temp2_I_sg, temp2i_sb, 
ray_top ); 

{ SET THE I t FIELD IN THE SOURCE RAY } 

temp2 I tr := o intensity_red; 
temp2 I tg := o intensity_green; 
temp2_I_tb := o_intensity_blue; 

{ RESTORE THE STACK } 

push( temp2_ray_type, 

temp2 ray_origin_x, temp2_ray_origin_y, temp2_ray_origin_z, 

temp2 ray vector x, temp2 ray vector y, temp2 ray vector _z, 

remp2 ray stvpe, 

temp2 jntersection rt ag, 

temp2_obj_idx, 

temp2 subobj idx, 

temp2_cpartidx, 

temp2_polygon_idx, 

temp2 intersection x, temp2 intersection y, temp2 intersection z, 
temp2_d, 

temp2_I_tr, temp2i_tg, temp2i_tb, 
temp2 I sr, temp2_I sg, temp2 I_sb, 
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ray_top 



); 



push( temp l_ray type, 

templ_ray_origin_x, templ_ray_origin_y, templ_ray_origin_z, 

tempi ray origin x, tempi ray origin y, tempi ray_origin_z, 

tempi _ray _stype, 

tempi _incerseccion_tlag, 

tempi objjidx, 

templ_subobj_idx, 

tempi c part _idx, 

tempi _polygon_idx, 

tempi intersection x, tempi intersection y, templ_intersection_z, 
templ_d, 

templ_I_tr, templ_I_tg, templ_Ijtb, 
templ_I_sr, templ_I_sg, templ_I_sb, 
ray_top ); 

end; 

end; 



{.PA} 

| **x**:«************x***********x** YlAll^ ************************ 

| ************************************************************************** 
begin 



*" } 
x«*xxxx***x i 



{ INPUT FILE } 
assign (sysin, ’pictures. pas’); 
reset (sysin); 



{ OUTPUT FILE } 
assign (outfile, ’picO.dta’); 
rewrite (outfile); 



{ 



************************************************************************** 



} 



{ SET UP COUNTERS } 
obj_cntr := 1; 
subobj_cntr := 1; 
light_cntr := 1; 
opart cntr := 1; 
poly cntr 1; 

vertice ^ncr := 1; 



{ CREATE POINTERS TO RECORDS } 
new(obj_curr); 
new(subobj_curr); 
new(cpart_curr); 
new(light jcurrent); 
new(polyjcurr); 
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{ *** READ IN DATA FILE *** } 



{ READ IN PICTURE RECORD } 
readln (sysin, picture. view_position_x); 
readln (sysin, picture. view_position_y); 
readln (sysin, picture. view_position_z); 
readln (sysin, picture. background_color_r); 
readln (sysin, picture.background_color_g); 
readln (sysin, picture, back ground _color_b); 
readln (sysin, picture.screen_max_x); 
readln (sysin, picture.screen_max_y); 
readln (sysin, picture. ambient_r); 
readln (sysin, picture. ambient_g); 
readln (sysin, picture. ambient b); 
readln (sysin, picture. nonzero); 
readln (sysin, picture. global refraction_index); 
readln (sysin, picture. num_lights); 
picture. lights := light jcurrent; 

while picture. num lights > 0 do begin 

{ READ IN LIGHT DATA } 

readln (sysin, light_current "(light_cntr].I_r); 
readln (sysin, light_current "(light_cntrj.I_g); 
readln (sysin, light_current "(light_cntr].I_b); 
readln (sysin, light_current "(light_cntr].light_x); 
readln (sysin, light current "(light cntr]. light y); 
readln (sysin, light current "(light cntr]. light z); 
readln (sysin, light_current "(light_cntr].dimensionl); 
readln (sysin, light_current "(light cntr].dimension2); 
light_cntr := light_cntr •+■ 1; 
picture. num_lights := picture. num_lights - 1; 
end; 

readln (sysin, picture. num_objs); 
object loop cnt := picture. num objs; 
picture. objects := obj curr; 



while object_loop_cnt > 0 do begin 

{ READ IN OBJECT DATA } 

readln (sysin, obj curr "(obj cntr]. opcode); 
readln (sysin, obj _curr "(obj cntr]. obj bsphere radius); 
readln (sysin, obj curr~ r obj cntri.obj bsphere x); 
readln (sysin, obj_curr "jobj _cntri.obj bsphere_y); 
readln (sysin, obj curr"|obj cntr). obj bsphere z); 
readln (sysin, obj_curr "(obj_cntr].num _sub_objects); 
subobj loop cnt := obj curr"[obj cntr]. num sub objects; 
picture. objects "(obj cntr]. sub objects := subobj curr; 

while subobj loop cnt > 0 do begin 

{ READ IN SUBOBJECT DATA } 

readln (sysin, subobj_curr "(subobj jcntr]. subobj _type); 
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readln (sysin, subobj_curr" subobj _cntrj.sub_bsphere_radius); 
readln (sysin, subobj jcurr " subobj_cntr].sub_bsphere_x); 
readln (sysin, s'ubobj curr" subobj cntrKsub bsphere y); 
readln (sysin, subobj_curr "[subobj cntrj.sub bsphere z); 
readln (sysin, subobj curr"subobj cntrj.num common parts); 
cpart loop cnt := subobj _curr "(subobj cntrl.num common parts; 
picture. objects " obj_cncr .suo_objects "Jsubobj_cncri. common parts : = 
cpart curr: 

while cpart_Joop_cnt > 0 do begin 
{ READ IN COMMON PART DATA } 

readlnisysin, cpartjrurr " opart cntr’.K ar); 
readln (sysin, cpart curr " cpart cntrj.K ag); 
readln (sysin, cpart j:urr "cpart jcntrj.K_ab); 
readln(sysin, cpart__curr " cpart^cntrj.K _dr); 
readln(sysin, cpart_curr " cpart_cntr].K_dg); 
readln(sysin, cpart_curr "[cpart cntrj.K db); 
readln(sysin, cpart curr " cpart cntrj.K sr); 
readlnisysin, cpartjcurr" cpart_cntrj.K_sg); 
readln(sysin, cpart_curr "cpart_cntrj.K_sb); 
readln (sysin, cpart_curr "cpart cntrj.Kjtr); 
readln (sysin, cpart_curr " cpart__cntr].K_tg); 
readln(sysin, cpart_curr" cpart _cntr!.K_tb); 
readlnisysin, cpart_curr " cpart cntrj.obj refraction__index); 
readln (sysin, cpart curr " cpart_cntr].obj_phong_exp); 

{ CHECK TO SEE IF SUBOBJECT IS A SPHERE OR A POLYGONAL OBJECT } 
if (subobj _curr"[subobj_cntr].subobj_type = 1) then begin 

readln(sysin, cpart_curr "(cpart_cntr].num_polygons); 
poly_loop_cnt := cpart_curr "[cpart _cntr].num_polygons; 
pict ure. objects "[obj_cntr].sub_objects "[subobj_cntr]. 
common_parts "[cpart_cntrj. polygons := poly_curr; 

while poly loop cnt > 0 do begin 
{ READ IN POLYGON DATA } 

readln(sysin, poly_curr "(poly_cntr].num_vertices); 
vertice loop cnt := poly curr "(poly cntrj.num vertices; 

while vertice loop_cnt > 0 do begin 
{ READ IN VERTICE DATA"} 

readln(sysin. poly curr " f poly cntrj. vertice xfvertice cntr]); 
readlnisysin. poiv curr"'poiy cntri. vertice yfvertice cntr;); 
reaaln|svsin. poly curr "'poly cntr'.vprtice zivertice cntr); 
vertice_cntr := vertice_cntr -r i; 
vertice loop cnt := vertice_Joop cnt - 1; 
end; 

vertice_cntr := 1; 

readln(sysin, poly_curr "[poly_cntr].surface_normal_x); 
readln(sysin, poly curr "(poly_cntr].surface_normal_y); 
readln(sysin, poly curr "[poly cntr]. surf ace normal z); 
poly_cntr := poly_cntr + 1; 
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poly_loop_cnt : — poly_loop_cnt - 1 
end; 

poly_cntr := 1; 
new(poly_curr); 

end; { IF } 

cparfc cntr := cpart_cntr 4- 1; 
cpart loop cnt cpart loop_cnt - 1 

end; 

cpart _cntr 1; 

new(cpart curr); 

subobj _cntr ;= subobj_cntr + 1; 

subobj loop cnt := subobj loop cnt - 1 

end; 

new (subobj _curr); 
subobj cntr ;= i; 
obj_cntr := objjcntr + 1; 
object_loop_cnt := object _loop_cnt - 1 



end; 

{.PA} 

{ SET RAY STACK POINTER } 

ray top := nil; 

pixel_z := initial_pixel_z; 

{ RASTER SCAN LOOP } 

for pixel_y := 1 to pict ure. screen _max_y do begin 

for pixel x := 1 to picture. screen max x do begin 

{ DETERMINE VIEW RAY DIRECTION } 
x := pixel_x - picture.view_position_x; 
y pixel y - picture. view position y; 
z := pixel_z - picture. view_position_z; 
dist := sqrtfsqr(x) 4- sar(y) 4* sqr(z)); 

{ CONVERT IT TO A UNIT VECTOR } 
umtx -= x dist; 
unity := y / dist; 
unitz := z / dist; 

{ INITIALIZE VIEW RAY } 

initial ray type := view; 
initial ray origin_x := pixel^x; 
initial^ray_origin_y := pixel y; 
initial ray origin z := pixel z; 
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initial_ray_vector_x := unitx; 
initial ray vector y := unity; 
initial_ray_vector_z := unitz; 
initial_ray_stype := none; 
initial_intersection_flag := false; 
initial objidx := 0; 

initial subobj idx := 0; 

initial_cpart_idx := 0; 

initial_polygon_idx := 0; 

initial intersection x := 0.0; 

initial_intersection_y := 0.0; 

initial_intersection_z := 0.0; 

initial_d := 0.0; 
initial_I_tr := 0.0; 
initial_I_tg := 0.0; 
initial_I_tb := 0.0; 
initial_I_sr := 0.0; 
initial_I_sg := 0.0; 
initial I_sb := 0.0: 

ray generation_number := 0; 

push( initial ray_type, 

initial ray origin_x,initial_ray_ongin_y,initial_ray_origin jz, 

initial_ray_vector_x,initial_ray_vector_y,initial_ray_vector_z, 

initial ray stype, 

initial intersection_flag, 

initial obj_idx, 

initial_subobj_idx, 

initial_cpart_idx, 

initial_polygon_idx, 

initial intersection x, 

initial_intersection_y, 

initial intersection z, 

initial_d, 

initial I tr, initial_I_tg, initia.1 I tb, 

initial I sr, initial I sg, initial I sb, 

ray_top ); 



{ BEGIN RAY TRACING LOOP } 

reDeat 

pop current ray type. 

current ray origin x, current ray_origin y, current ray origin z, 

current _ray_vector_x,current_ray _yector_y,current_ray_vector_z, 

current_ray_stype, 

current intersection flag, 

current _obj_idx, 

c urre n t _su bobj _id x , 

current cpart idx, 

current_polygon_idx, 
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current_intersection_x, 
current_intersection_y, 
current_intersection z, 
current d, 

current I tr, current_I_tg, current I tb, 
current I sr, current I sg, current I sb, 
ray^top ); 



if (current intersection flag) then begin 

{ 

IF THIS FLAG IS SET THEN THIS RAY HAS ALREADY BEEN THROUGH THE RAY 
TRACING PROCESS AND HAS HIT AN OBJECT. 



} 



cpart_path := picture. 

objects "[current obj idx]. 

sub objects "[current__subobj idx]. 

common parts; 



subobj path := picture. 

objects "[current_obj_idxj. 
sub_objects; 

if subobj_path " [current_subobj_jdx|. subobj _type = 0 then begin 

{ 

IF THE SUBOBJECT TYPE IS A SPHERE THEN THE SURFACE NORMAL AT THE POINT 
OF INTERSECTION MUST BE CALCULATED SINCE IT CAN NOT BE STORED. THE SURFACE 
NORMAL IS DETERMINED FOR THE PLANE TANGENT TO THE SPHERE AT THE 
INTERSECTION POINT. 

} 

surface_normal_x := current_intersection_x - 

subobj _path " [current_subobj idx] . 

sub _bsphere_x; 

surface_normal_y := current_intersection_y - 

subobj_path "[current subobj idx]. 
sub_bsphere y; 

surface_normal_z := current_intersection_z - 

subobj path "[current subobj idx]. 
sub_bsphere z; 

dist := f sqrtf sqrfsurface normal x) +- 
sqr( surface normal y) -+- 
sqrfsurface _normai z ) ) ) ; 

{ THIS RAY IS THEN CONVERTED INTO A UNIT VECTOR } 
surf ace _normal_x := surface_normal_x / dist; 
surface_normal y := surface_normal_y / dist; 
surface normal z := surface normal z / dist; 
end 

else begin 

{ 
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IF THE SUBOBJECT IS A POLYGONAL OBJECT THEN THE SURFACE NORMALS FOR EACH 
OF THE POLYGONS OF WHICH IT IS COMPOSED IS RETRIEVED FROM IT’S RECORD. 

} 

surface_normal_x := cpart_path "|current_cpart_idx]. 

polygons "(current polygon idx]. 
surface normal x: 

surface_normal_y := cpart_path "|current_cpart_idxj. 

polygons "[current polygon_idx]. 
surface normal y; 

surface normal z cpart path "[current cpart idx]. 

polygons "[current polygon idx]. 
surface_normal_z; 

end; 

{ 

PROCEDURE FOR DETERMINING THE INTENSITY OF LIGHT AT EACH INTERSECTION 
POINT. 

} 

12 (current_d, 

current_I_tr, current_I_tg, current_I_tb, 

current^! _sr, current __I s g , current _I_tb, 

picture, num flights, 

picture. ambient_r. picture. ambient jg, picture. ambient_b, 

cpart_path "[current cpart idx .K ar, 

cpart_path "[current_cpart_idx].K_ag, 

cpart path "[current cpart idxJ.K ab, 

cpart jpath "[current_cpart_idx\K_dr, 

cpart_path " [current _cpart_idx].K_dg, 

cpart_path "[current cpart idxj.K db, 

cpart path "[current cpart idxj.K sr, 

cpart __path "(current_cpart_idx].K_sg, 

cpart _path "[current_cpart_idx].K_sb, 

cpart _path "[current_cpart_idx].K_tr, 

cpart_path "(current _cpart_idx].K_tg, 

cpart _path "[current_cpart_idx].K_tb, 

current_ray type, 

ray_top, 

picture. lights, 

picture. objects, 

cpart path "-current cpart idx .obj phong exp, 

current rav vector x. 

current rav vector y, 

current ray vector z, 

current intersection x, 

current_intersection_y, 

current_intersection_z, 

surface_normal_x, surface normal y, surface normal z, 

ray_generation_number, 

intensity red, 

intensity green, 
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intensity blue); 



end 

else begin 



{ CHECK FOR POSSIBLE INTERSECTIONS OF CURRENT RAY WITH OBJECTS IN SCENE } 

check for intersection(current ray vector x, 
current _ray_vector_y, 
current_ray vector_z, 
current ray origin x, 
curren t_ray_origin_y, 
current ray origin z, 
picture. objects, 
current_intersection_x, 
current intersection y, 
current intersection z, 
curren t _obj _idx, 
current_subobj_idx, 
current_cpart idx, 
current polygon idx, 
current_intersection_flag); 



{ SET UP PATHNAMES TO USE AS SHORTHAND } 

subobj_path := picture. 

objects "[current_obj_idx]. 
sub objects; 

cpart_path := picture. 

objects "[current obj idx). 
sub_objects "[current _subobj_idx]. 
common_parts; 

if subobj_path "(current jsubobj_idx].subobj_type = 0 then begin 
surface_normal_x := current_intersection_x - 

subobj path "(current subobj idx]. 
sub bsphere x; 

surface_normal_y := current_intersection_y - 

subobj _path " [curren t_subobj _idx] . 
sub bsphere y; 

surface normal z := current intersection z - 

subobj path " [curren t subobj idx]. 

sub bsphere z\ 

disc ( sqrt( sqr(surface normal x) - 
sqr(surface normal yj -r 
sqr(surface_normal_z))); 

surface normal_x := surface_normal_x / dist; 
surface normal y := surface normal y / dist; 
surface normal z := surface normal z / dist; 
end 

else begin 
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surface_normal_x := cpart_path ~[current_cpart_idx]. 

polygons "[current polygon idx] . 
s urf ac e _n orm al x ; 

surf ace _normal_y := cpart_path "[current_cpart Jidx]. 

polygons "fcurrent polygon _idxl. 
surface normal y; 

surface normal_z := cpart_path "[current_cpart_idx]. 

polygons "[current_polygon_idx]. 
surface normal z; 

end; 

if (current intersection_fiag) then begin 

{ IF THERE HAS BEEN AN INTERSECTION THEN CONTINUE TRACING THE RAY } 

if (stack exceeded ( ray generation number, 

maximum size jof_stack)) then begin 

{ IF THE STACK IS ALREADY FULL THEN CALCULATE INTENSITY AT LAST NODE } 
16 (current_d, 

current I_tr, current I tg, current I tb, 
current I sr, current I sg, current I tb, 
picture. num_lights, . 
picture. ambient_r, 
picture. ambient_g, 
picture. ambient_b, 

cpart path "[current cpart idx].K ar, 
cpart_path "[current_cpart idx].K_ag, 
cpart_path " [current _cpart_idx].K_ab, 
cpart _path "[current_cpart_idx] .K_dr, 
cpart_path "[current cpart idx].K dg, 
cpart_path "[current_cpart_idx].K_db, 
cpart path "[current cpart idxJ.K sr, 
cpart_path " [current_cpart_idx] .K_sg, 
cpart_path "[current_cpart_idx].K_sb, 
cpart path "[current cpart idx].K tr, 
cpart_path "[current_cpart_idx].K_tg, 
cpart_path "[current_cpart_idx].K_tb, 
current _ray_type, 
ray_top, 
picture. lights, 
picture. objects. 

cpart path "'current __cpart idxi.obj phong exp, 

current_ray_vector_x, 
current_ray_vector_y, 
current ray_vector z, 
current intersection x, 
current_intersection_y, 
current_intersection z, 
surface normal x, 
surface_normal_y, 
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surface_normal_z, 
intensity red, 
intensity_green, 
intensity_blue) 

end 

else begin 

{ 

IF THERE WAS AN INTERSECTION AND THE STACK WAS NOT FULL THEN CALCULATE 
THE DISTANCE BETWEEN THE RAY’S ORIGIN AND POINT OF INTERSECTION AND PLACE 
THE RAY BACK ON THE STACK. 

} 



current_d := (sqrt(sqr(current_intersection_x - 
current_ray_origin_x)) H- 
(sqr(current intersection_x - 
current_ray_origin_x)) -f 
(sqr(current intersection x - 
current _ray origin x) ) ) ; 

push( current_ray_type, 
current_ray_origin_x, 
current ray origin y, 
current _ray_origin_z, 
current ray vector x, 
current ray vector y, 
current _ray_vector_z, 
current _r ay _stype, 
current_intersection_flag, 
current obj idx, 
current_subobj_idx, 
current cpart idx, 
current polygon idx, 
current_intersection_x, 
current_intersection_y, 
current_intersection_z, 
current d, 
current_I_tr, 
current I tg, 
current_I_tb, 
current_I_sr, 
current_I_tg, 
current I fcb, 
ray_top ); 

{ DETERMINE IF A REFLECTED RAY WAS CREATED AND IF SO CALCULATE IT } 
calculate_reflected_ray (current _ray_vector_x, 
current_ray^vector__y, 
current ray vector z, 
surface normal x, 
surface_normal_y, 
surface normal z, 
reflected ray x, 
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reflected ray v, 
reflected ray z , 
reflected ray ); 



{ DETERMINE IF A REFRACTED RAY WAS CREATED AND IF SO CALCULATE IT } 
calculate refracted _ray (current ray v r ector x, 
current ray vector y, 
current ray vector z, 
surface normal x, 
surface_normal_y , 
surface_normal_z, 
cpart path "[current cpart idx]. 
obj_refraction_index, 
picture. global refraction index, 
refracted ray x, 
refract ed_ray_y, 
refracted ray z, 
refracted ray ); 

ray generation number := ray generation number + 1; 
if reflected ray then begin 



{ 

IF A REFLECTED RAY WAS CREATED THEN INITIALIZE IT AND PUSH IT ON THE 
STACK. 

} 

source ray type := current ray type; 

dist := ( sqrt(sqr(reflected_ray_x) + 
sqr(reflected_ray_y) + 
sqr(reflected_ray _z) ) ) ; 

{ CONVERT REFLECTED RAY TO A UNIT VECTOR } 

unitx := reflected ray x / dist; 
unity := reflected_ray_y / dist; 
unitz := reflected ray z / dist; 

initial_ray_type := reflected; 
initial ray origin x := current intersection x; 

initial ray ^origin y currentintersection y; 

initial _ray_origm z current ^intersection z: 

initial rav vector x :*= unitx; 
mitiai_ray_vector_y := unity; 
initial_ray_vector_z := unitz; 
initial_ray_stype := source_ray^type; 
initial_intersection_flag := false; 
initial_obj idx := 0; 

initial_subobj_idx := 0; 

initial^cpart^idx := 0; 

initial_polygon idx := 0; 
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initial intersect ion _x := 0.0 

initial_intersection_y := 0.0 

initial intersection_z := 0.0 

initial_d := 0.0; 

initial_I_tr := 0.0; 

initial I_tg := 0.0; 

initial [ tb := 0.0; 

initial ^1 sr := 0.0; 

initialling := 0.0; 

initial I sb := 0.0; 



push( initial_r a y_ fc yp e i 
initial_r a y J° r igi n _ x ) 
initial _ray __origin_y» 
initial_r a y_ or igi n _ z ? 
initial ray_v ector _ x > 
initial ray vector_y 5 
initialj* a y _v ec ^ or _ z > 

initial ray stype, 
initial intersection flag, 
initial_obj_id x , 
initial _subobj_idx, 
initial _cp ai *t_idx, 
initial_p°lyg° n _i ( ix, 

initial_intersection_Xj 

initial_intersection_y> 

initial intersection z, 
initial d, 

initial I_tr, initial_I_tgj initial_I tb, 

initial I sr, initial I_sg> initial l_sb, 
ray_top ); 



end; 



{ 



if refracted ray then begin 



IF A REFRACTED RAY WAS CREATED THEN INITIALIZE IT AND PUSH IT ON THE 
STACK. 



} 



source ray type current ray type; 
dist ( sqrt(sqr(refracted ray x) -+■ 

sqr(refracted_r a y_y) + 

sqr(refracted_r a y z))); 

{ CONVERT IT TO A UNIT VECTOR } 

unitx := refracted _r a y_ x / <dist; 
unity := refracted ray y / dist; 
unitz := refracted_r a y_ z / dist; 
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imtial_ray_type := refracted; 
initial_ray_origin x := current intersection x; 
initial_ray_origin_y := current_intersection_y; 
initial ray_origin z := current intersection z; 



initial_ray_vector_x := unitx; 
initial_ray_vector y := unity; 
initiai_ray_veccor__z unitz; 

initial_ray_stype := source_ray_type; 
initial intersection flag := false; 
initial_obj_idx := 0;- 

initial_subobj_idx := 0; 

initial _cpart_idx := 0; 

initial_polygon_idx := 0; 

initial_intersection_x := 0.0; 

initial_intersection_y := 0.0; 

initial_intersection_z := 0.0; 

initial_d := 0.0; 
initial _I_tr := 0.0; 

initial_I_tg := 0.0; 

iniual_I_tb := 0.0; 

initial_I_sr := 0.0; 

initial_I_sg := 0.0; 

initial I sb := 0.0; 



push( initial ray type, 
initial_ray_origin x, 
initial_ray_origin_y, 
initial_ray origin z, 
initial_ray_vector_x, 
initial_ray vector y, 
initial ray vector z, 
initial_ray_stype, 
initial intersection flag, 
initial_obj_idx, 
initial_subobj_idx, 
initial cpart idx, 
initial_polygon_idx, 
initial_intersection_x, 
initial_intersection_y, 
initial_intersection_z, 
initial d, 

initial I tr. initial I tg, initial I tb, 
initial Ijr. initial _l_sg, initial ^[_sb, 
ray_top j; 



end; 

end; 

end 

else begin 
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if (current_ray_type = view) then begin 

{ 

IF THERE WAS NO INTERSECTION AND THE CURRENT RAY IS THE VIEW RAY THEN SET 
THE OUTPUT INTENSITY TO THE BACKGROUND INTENSITY. 

} 

intensity red := picture. background color r; 
intensity_green := picture. oackground color g; 
intensity_blue := picture. background_color b; 

end; 

end; 

end; 

until (stack_empty(ray_top)); 

{ OUTPUT THE FINAL INTENSITY } 
display_plxel( intensity_red, 
intensity_green, 
intensity Jblue, 
pLxel_x, 
pixel_y, 
pixel_z ); 

end 

end; 

close (sysin); 
close (output); 

end. { MAIN } 
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APPENDIX B - INPUT FILE 



100 x/ view position { PICTURE } 

100 /y/ 

1000 /z/ 

0.0 /red/ background light 
0.0 / green/ 

1.0 /blue/ 

200 /x/ screen 

200 /y/ 

1.0 /red/ ambient intensity 

1.0 /green/ 

1.0 /blue/ 

1.0 /global refraction index/ 

1 numlighcs/ 

1.0 /red/ intensity of light source 

1.0 /green/ 

1.0 /blue/ 

0.0 /x/ position of light source 

20.0 /y / 

0.0 jzj 

0.0 /dimensionl/ 

0.0 /dimension2/ 

3 /num_objects/ 

9999 /opcode/ { OBJECT 1 } 

35.0 /radius of object's bounding sphere/ 

110.0 /x/ center of bounding sphere 

0.0 /y/ 

-50.0 /z/ 

1 /number of subobjects/ 

1 /subobject type/ {SUBOBJECT 1} 

35.0 /radius of subobjects bounding sphere/ 

110.0 /x/ center of bounding sphere 

0.0 /y/ 

-50.0 /z / 

1 /num-common-parts/ { COMMON PART 1 } 



0.8 


/Ka-red/ 


ambient coefficient 


0.0 


/Ka- green/ 




0.0 


Ka-biue/ 




0.3 


Kd-red/ 


diffuse coefficient 


0.0 


i Kd-green/ 




0.0 


/Kd-blue/ 




0.8 


/Ks-red / 


specular coefficient 


0.8 


/Ks- green/ 




0.8 


/Ks-blue/ 




0.0 


/Kt-red/ 


transmission coefficient 


0.0 


/Kt-green/ 




0.0 


/Kt-blue/ 
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0.0 

200 

6 

4 

90.0 

20.0 
- 20.0 

90.0 
- 20.0 
-30.0 

130.0 
- 20.0 
-30.0 

130.0 

20.0 
-30.0 
0.0 
0.0 
1.0 

4 

90.0 

- 20.0 

-30.0 

90.0 
- 20.0 
-70.0 

130.0 
- 20.0 
-70.0 
130.0 
- 20.0 
-70.0 
0.0 
-1.0 
0.0 

4 

130.0 

20.0 
-70.0 

130.0 
- 20.0 
-70.0 

90.0 
- 20.0 
-70.0 

90.0 

20.0 
-70.0 
0.0 
0.0 
-1.0 
4 



/obj-refraction-index/ 

/obj-phong-specular-exponent/ 

/num-polygons/ 

/num-vertices/ { POLYGON 1 } 
/poly 1 ptl/ 



/polyl pt2/ 



/polyl pt3/ 



/polyl pt4/ 



/polyl surface normal/ 



/num-vertices/ { POLYOGN 2 } 
/poly2 ptl/ 



/poly 2 pt2/ 

/poly2 pt3/ 

/poly2 pt4/ 

/poly2 surface normal/ 

/num-vertices/ {POLYGON 3} 
/poly3 ptl/ 

/poly3 pt2 / 



poly 3 pr,3/ 



/poly3 pt4/ 



/poly3 surface normal/ 
/num-vertices/ { POLYGON 4 } 



130.0 /poly4 pt 1/ 

20.0 

-30.0 

130.0 /poly4 pt2/ 

20.0 

-70.0 

90.0 /poly4pt3/ 

20.0 

-70.0 

90.0 /poly4 pt4/ 

20.0 

-30.0 

0.0 /poly4 surface normal/ 

1.0 

0.0 

4 /num-vertices/ { POLYGON 5 } 

130.0 /poly 5 ptl/ 

20.0 
-30.0 

130.0 /polyS pt2/ 

- 20.0 

-30.0 

130.0 /poly5 pt3/ 

- 20.0 

-70.0 

130.0 /poly5 pt4/ 

20.0 

-70.0 

1.0 /poly 5 surface normal/ 

0.0 

0.0 

4 /num-vertices/ { POLYGON 6 } 

90.0 /poly6 ptl/ 

20.0 

-70.0 

90.0 /poly6 pt2/ 

- 20.0 

-70.0 

90.0 /poly6pt3/ 

- 20.0 

-30.0 

90.0 /polv6 pt4/ 

20.0 
-30.0 

-1.0 /poly6 surface normal/ 

0.0 

0.0 

9999 /opcode/ *** OBJECT 2 *** 

175 /radius of the objects bounding sphere/ 

100.0 /x/ center point of the bounding sphere 

-100.0 /y/ 

-100.0 /z/ 
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1 /num-subobjects/ 

1 /subobject-type/ { SUBOBJECT 1 } 

175 /radius of subobjects bounding sphere/ 

100.0 /x/ center of bounding sphere 
-100.0 /y / 

-100.0 h/ 

1 /number common parts/ { COMMON PART 1 } 

0.0 /red/ Ka ambient coefficient 
0.7 /green/ 

0.0 /blue/ 

0.0 /red/ Kd diffuse coefficient 
0.7 /green/ 

0.0 /blue 

0.8 /red/ Ks specular coefficient 
0.8 /green/ 

0.8 /blue/ 

0.0 /red/ Kt transmission coefficient 
0.0 /green/ 

0.0 /blue/ 

0.0 /objects refraction index/ 

200 /Phong’s specular exponent/ 

1 /number of polygons/ 

4 /number of vertices/ { POLYGON 1 } 

0.0 /polyl pt 1 / 

20.0 
- 200.0 

0.0 /polyl pt 2/ 

0.0 

0.0 

200.0 /polyl pt3 / 

0.0 

0.0 

200.0 /polyl pt4/ 

20.0 
- 200.0 

0.0 /polyl surface normal/ 

0.99 

0.1 

9999 /opcode/ * * * * OBJECT 3 * * * * 

40 /radius of objects bounding sphere/ 

140 /x/ center of bounding sphere 

30 /y / 

-150 Izj 

1 /number of subobjects/ 

0 /subobject type/ 

40 /radius of subobjects bounding sphere/ 

140 /x/ center of bounding sphere 

30 /y/ 

-150 /z/ 

1 /number of common parts / { COMMON PART 1 } 
0.5 /red/ Ka ambient coefficient 

0.0 /green/ 
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diffuse coefficient 



0.5 /blue/ 

0.5 /red/ Kd 
0.0 /green/ 

0.5 /blue/ 

0.8 / red/ Ks specular coefficient 

0.8 1 green/ 

0.8 blue/ 

0.0 /red/ Kt transmission coefficient 

0.0 /green/ 

0.0 /blue/ 

0.0 /refraction index for object/ 

200 / Phong ? s specular exponent/ 
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